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ABSTRACT 

The  two-phase  flow  in  a convergent-divergent  nozzle  was  analyzed. 
The  model  used  is  that  of  one-dimensional  equilibrium  conditions, 
uniform  temperature  in  a cross-section  and  zero  wall  friction. 

The  liquid  phase  is  assumed  to  consist  of  uniformly  distributed 
spherical  droplets  of  equal  diameter,  whose  size  is  determined 
by  the  limiting  Weber  number  breakup  criterion.  Smaller  droplets 
may  be  fed  to  the  nozzle  initially.  Formulas  for  mass-flux 
density,  thrust  and  overall  efficiency  are  given.  The  results 
provided  by  a computer  program  include  the  slip  ratio  (i.e. 
velocity  ratio  of  the  two  phases)  along  the  nozzle  axis  and 
the  corresponding  variation  of  the  droplet  size  and  the  local 
friction  losses  due  to  droplet  drag.  The  cone  angle  of  the 
nozzle  may  be  prescribed  or  it  may  be  determined  from  a selected 
slip  ratio  gradient. 

Initial  conditions  of  slip  ratio  and  droplet  size  required  were 
investigated  in  detail.  The  ratio  of  droplet  size  to  nozzle 
diameter  is  found  to  vary  with  nozzle  size  because  of  Reynolds 
and  Weber  number  effects. 

The  slip  model  is  found  to  predict  a much  higher  mass-flux 
density  at  the  throat  than  that  given  by  the  I.H.E.  (Isentropic, 
Homogeneous,  Equilibrium)  model.  Maximum  mass-flux  densities 
and  minimum  droplet  sizes  at  the  throat  are  outputs. 

The  two-dimensional  two-phase  flow  in  impulse  turbine  buckets  was 
also  investigated.  The  parameters  that  control  the  deviation 
angle  of  the  droplets  from  the  steam  or  gas  path  are  identified. 
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INTRODUCTION 

Two-phase  nozzles  are  employed  to  accelerate  a cloud  of  liquid 
droplets  suspended  in  a gaseous  carrier  to  impart  kinetic  energy 
to  the  liquid  phase  for  the  purpose  of  driving  a two-phase  or  liquid 
turbine.  Large  available  pressure  heads  result  in  relatively 
modest  discharge  velocities  of  the  mixtures  because  the  droplets 
effectively  increase  the  molecular  weight  of  the  mixture.  Since 
the  sonic  velocity  in  the  mixture  is  also  lowered,  the  spouting 
velocity  of  the  nozzle  is  generally  supersonic,  resulting  in  a 
converging-diverging  geometry. 


The  performance  of  a two-phase  nozzle  is  expected  to  increase 
with  decreasing  droplet  size.  A mathematical  model  has  been 
developed  to  study  the  quantitative  relationship  between  nozzle 
design,  droplet  size,  and  relative  velocity  or  slip  between  the 
liquid  and  gaseous  phases  and  the  effects  of  irreversibilities. 
That  is  an  important  consideration,  since  the  characteristics  of 
the  mixture  will  subsequently  affect  the  turbine  design  and 
performance . 

The  simplest  model  would  be  the  so-called  Isentropic  Homogeneous 
Equilibrium  (IHE)  model1  which  does  not  include  slip  and  therefore 
friction  losses  or  heat-transfer  across  the  nozzle  wall,  nor  does 
it  admit  any  metastable  states.2’3 


\ 


The  problems  of  supersonic  nozzle  design  are  the  determination  of 
the  following: 

1.  Critical  mass-flux  density  or  the  size  of  the  throat 
area  for  a given  mass  flow  rate. 

2.  The  required  area  ratio  (discharge  area/throat  area). 


3.  Predicted  performance  (efficiency)  as  a function  of 
nozzle  geometry,  initial  pressure,  initial  vapor 
fraction,  approach  velocity  (and  possibly  time  available 
for  droplet  breakup) . 
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The  first  two  may  be  well  approximated  for  higher  vapor  fractions 
(i.e.,  x > 0.2)  by  the  I HE  model.  However,  the  region  of  most 
interest  in  two-phase  turbine  applications  is  that  for  which  the 
IHE  model  gives  the  poorest  predictions.  Thus  a more  sophisti- 
cated model  that  considers  slip  must  be  developed  to  enable 
nozzle  predictions  to  be  made  for  two-phase  turbines. 

The  complexity  of  such  an  analysis  depends  on  whether  or  not 
further  evaporation  occurs  during  the  course  of  the  expansion, 
which  changes  the  vapor  fraction. 

Case  I,  with  constant  vapor  fraction,  is  approached  by  two  compo- 
nent systems  consisting  of  a mixture  of  initially  saturated  or 
superheated  steam  and  an  oil  of  low  vapor  pressure. 

Also,  mixtures  of  air  with  water,  initially  at  room  temperature 
and  elevated  pressure,  fall  into  that  category. 


Case  II,  with  evaporation,  is  represented  by  single  component 
systems  exemplified  by  expanding  wet  steam. 

Case  III,  which  involves  evaporation  in  a two-component  system, 
is  found  for  example  in  an  initially  pressurized,  heated  mixture 
of  two  liquids,  to  be  expanded  such  that  one  liquid  flashes  into 
vapor.  That  case  is  not  treated  here  in  detail,  since  the 
emphasis  was  placed  on  Case  II,  the  analysis  of  wet  steam  nozzle 
flow  with  droplets  trailing  the  steam.  The  application  of  the 
analysis  to  Case  I with  constant  vapor  fraction  which  is  realized 
for  example  in  a mixture  of  water  droplets  in  air  is  relatively 

easily  made  (it  is  a special,  simpler  case)  but  is  not  pursued 
further  at  this  time. 

Biphase  Energy  Systems  presently  has  a computer  code  in  its 
library  that  can  predict  nozzle  performance  in  the  region  of 
interest.  However,  the  code  requires  the  input  of  a pressure 
profile  along  the  nozzle  axis,  with  the  nozzle  geometry  being 


an  output.  This  restriction  makes  the  off-design  study  of  nozzle  8 

performance  cumbersome.  In  particular,  it  is  difficult  to  study 
the  performance  of  straight  conical  nozzles  (which  are  easy  to 
make)  with  the  present  code.  In  addition,  the  existing  program 
includes  heat- transfer  effects  between  the  two  phases  which  makes 
the  analysis  more  cumbersome.  Independent  analysis  and  experimen- 
tal results  show  that  for  small  droplets  the  heat  flux  is  very 
high  and  therefore  a more  simplified  approach  may  concentrate  on 
the  more  important  influences  on  the  droplet  formation  in  a two- 
phase  nozzle  with  slip. 

Such  a simplified  analysis  of  the  physical  phenomena  occurring  in 
two-phase  flow  is  presented  in  Part  I. 

While  the  progression  of  droplet  breakup,  of  local  friction  losses 
between  the  slower  moving  droplets  and  the  driving  steam  and  the 
variation  of  the  slip  (or  velocity)  ratio  is  studied  by  means  of 
the  computer  program  which  integrates  numerically  the  differential 
equation  for  the  slip  ratio  step  by  step  along  the  expansion  path, 
certain  conclusions  can  be  directly  developed  from  the  equations 
about  the  conditions  at  the  nozzle  inlet  and  at  the  throat.  Some 
such  results  are  the  initial  slip  ratio  and  the  ratio  of  droplet 

$ 

size  to  nozzle  diameter  at  the  inlet  as  a function  of  nozzle  inlet 
size . 


At  the  throat,  upper  limits  of  the  mass-flux  density  can  be  deter- 
mined as  a function  of  temperature  and  slip  ratio  at  the  throat. 
Such  results  are  of  prime  interest  because  of  the  considerable 
under -prediction  of  the  mass-flow  rate  of  wet  steam  nozzles  with 
small  initial  vapor  mass  fractions  by  the  simpler  theories. 


It  will  be  seen  that  the  theory  given  here  indeed  predicts  higher 
mass-flux  density,  and  also  shows  that  small  nozzles  develop  smaller 
droplets  at  the  inlet,  although  somewhat  larger  than  would  be  given 
by  a direct  proportion  to  the  opening  diameter.  The  initial  slip 
ratio  is  correspondingly  larger  for  the  small  nozzle.  Lower  inlet 
temperature  and  vapor  fraction  also  make  the  achievement  of  good 
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nozzle  performance  a challenge.  The  study  of  the  conditions  at 
the  nozzle  inlet  will  determine  the  initial  droplet  size  that 
needs  to  be  fed  to  the  nozzle  in  order  to  achieve  certain  per- 
formance objectives.  The  question  of  how  the  preatomization 
could  most  effectively  be  achieved  is  not  treated  here.  A sep- 
arate bibliography  compiled  in  Appendix  I gives  an  idea  of  the 
work  already  done  in  that  area  as  well  as  other  topics  related 
to  two-phase  flow  in  turbines  like  erosion  characteristics. 

The  results  of  Part  I should  primarily  be  considered  as  work- 
ing tools  for  a design  study  of  a two-phase  engine.  The  formu- 
lations provided  allow  the  calculation  of  all  that  is  needed  in 
a basic  model  study  when  the  initial  vapor  fraction  and  the  slip 
ratio  are  given.  When  these  are  unknown,  they  may  be  found  by 
means  of  the  computer  program. 

In  Part  II,  two-phase  flow  in  the  turbine  is  considered.  A result 
of  the  nozzle  study  is  the  final  droplet  size  emerging  from  the 
nozzle.  That  leads  to  a consideration  of  whether  more  conven- 
tional turbines  may  be  used  effectively,  that  is  without  undue 
separation  and  impingement  of  droplets  on  the  bucket  walls.  The 
formation  of  a thin  liquid  film  on  the  buckets  is  associated  with 
considerable  friction  losses  because  of  the  unfavorably  small 
hydraulic  diameter  of  such  a film.  Special  bucket  shapes  that 
concentrate  the  film  to  a larger  hydraulic  diameter  may  alleviate 
the  problem  and  were  proposed  in  Reference  8 . The  deviation  of 
the  droplet,  path  from  the  path  of  the  gas  is  treated  in  Part  II. 

Besides  the  study  of  particle  trajectories,  a turbine  design 
would  have  to  be  adapted  to  accommodate  supersonic  two-phase  flow 
through  the  buckets  at  starting  and  design  conditions.  The 
reader  is  referred  to  items  8 to  13  in  the  Bibliography  (page 
T-41 , 42)  of  Reference  4. 
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The  results  of  Part  II  may  not  enable  final  conclusions  to  be 
drawn  concerning  the  general  feasibility  of  two-phase  impulse 
turbines.  However,  the  tools  given  will  form  the  basis  for  a 
specific  study  to  be  proposed. 

APPROACH  AND  RESULTS 

The  present  work  uses  first  principles  to  develop  the  basic 
relationships  that  influence  nozzle  performance.  The  analysis 
is  similar  to  that  presented  in  References  5 and  6 which  is 
very  general  when  the  loss  mechanism  is  excluded.  The  loss 
formulation  of  References  5 and  6 used  an  impact  (billiard  ball) 
model  between  fluid  elements  moving  at  different  velocities. 

The  impact  model  is  replaced  in  this  work  by  the  model  employed 
by  Elliott  and  Weinberg7  which  is  based  on  the  orderly  movement 
of  a cloud  of  spherical  liquid  droplets  of  equal  diameter  uni- 
formly dispersed  in  the  gaseous  phase.  Other  assumptions  used  in 
the  analysis  are: 

1.  The  fluids  are  in  chemical  and  thermal  equilibrium. 

2.  Flow  is  one-dimensional  in  the  direction  of  the  nozzle 
axis . 

3.  No  friction  at  the  wall. 

4.  Droplet  breakup  is  controlled  by  a limiting  Weber 
number  criterion. 

Since  Reference  5 was  based  on  the  expansion  of  wet  steam,  the 
same  single  component  fluid  was  used  here,  even  though  the  con- 
sideration of  Case  II  with  continuously  evaporating  water  droplets 
(and  therefore  changing  vapor  fraction)  was  more  complicated 
analytically  than  the  study  of  Case  I with  a constant  vapor 
fraction.  Also,  the  development  of  an  algorithm  for  wet  steam 
properties  was  considerably  more  involved  than  the  use  of  ideal 
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gas  properties.  However,  the  potential  practical  usefulness  of 
such  a representation,  particularly  in  a computer  simulation, 
which  includes  viscosity  and  surface  tension,  was  determined  to 
justify  the  additional  complexity. 

The  analysis  is  presented  in  detail  in  Part  I.  The  first  step 

was  the  formulation  of  the  basic  ’ifferential  equation  containing 

the  derivative  of  the  slip  ratio  with  respect  to  pressure.  A 

loss  mechanism  was  not  included  in  the  early  formulation,  but, 

rather,  the  relative  loss  term  (1  - n ) was  employed  as  a para- 

N 

meter  in  the  analysis.  The  simplified  relationships  (Equation 
(58)  in  Part  I)  were  programmed  on  the  computer  and  parameter  sweeps 
were  run  to  determine  the  initial  slip  ratio  and  droplet  size 
without  preatomization. 

The  next  step  was  to  formulate  the  relative  loss  based  on  the 
frictional  drag  of  the  trailing  droplets.  The  initial  slip  ratio 
obtained  from  the  simplified  analysis  was  used  to  predict  the 
relative  loss  of  energy  in  the  first  increment  of  the  nozzle. 

The  resulting  governing  differential  equation  is  given  by 
Equation  (90),  which  is  an  ordinary  equation  of  the  first  order. 

As  pointed  out  in  the  introduction,  a separate  program  was  written 
for  the  determination  of  the  initial  slip  ratio  K to  be  used  at 
the  inlet  to  the  nozzle.  A criterion  was  found  in  the  form  of 
dK/dp  = 0 (Equation  (90))  which  enables  a limiting  droplet 
size  dL  to  be  determined  for  an  assumed  slip  ratio  K^.  Which 
dL  and  should  then  be  used  is  a question  of  the  droplet  size 
that  can  be  provided  at  the  inlet.  One  approach  is  to  use  the 
Weber  Number  criterion  for  the  droplet  breakup  which  in  effect 
gives  another  curve  dw  as  a function  of  K.  The  intersection  of 
the  d and  d.  lines  gives  the  solution  at  the  nozzle  inlet  if  no 
preatomization  is  provided  for.  Another  approach  is  to  prescribe 
the  maximum  droplet  size  at  the  nozzle  inlet  by  the  use  of  a 
separate  preatomizer.  The  prescribed  droplet  diameter  is  then 
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set  equal  to  dL;  the  dL  curve  pertaining  to  the  selected  inlet 
diameter  will  then  give  the  initial  slip  ratio  K. 

The  above  study  also  provides  a basis  for  a study  of  similarity 
relationships  as,  for  example,  the  ratio  of  droplet  to  nozzle 
diameter  at  the  nozzle  inlet  for  geometrically  similar  nozzles. 
Reference  8 describes  an  inve  tigation  of  similarity  relations 
which  was  based  upon  small  slip  ratios.  Conclusions  were  drawn 
as  to  the  effect  of  vapor  fraction  and  initial  pressure  levels. 

The  loss  formulation  (Equation  (82)  and  (83))  given  in  Part  I 
shows  that  the  ratio  of  droplet  diameter  to  nozzle  diameter  is 
a significant  non-dimensional  parameter.  It  seems  that  a simple 
similarity  rule  were  possible  in  terms  of  the  initial  droplet 
size/nozzle  diameter  ratio  if  it  were  not  for  the  presence  of 
the  Reynolds  number  effect  on  the  droplet  drag  coefficient.  Also, 
use  of  the  Weber  number  criterion  for  the  droplet  breakup  intro- 
duces effects  which  cause  deviations  from  a simple  similarity 
rule  based  on  an  initial  droplet  size/nozzle  diameter  ratio. 
Details  are  given  in  Part  I. 

Besides  the  special  study  of  the  initial  conditions  required  at 
the  nozzle  inlet,  the  general  theory  presented  in  Part  I can  also 
be  used  to  establish  limits  of  the  mass-flux  density  at  the  nozzle 
throat . Conclusions  regarding  possible  droplet  size  can  also  be 
drawn.  Provided  the  slip  ratio  K is  decreasing  in  the  diverging 
part  of  the  nozzle,  the  droplet  size  would  be  increasing  according 
to  the  Weber  number  criterion.  However,  it  is  reasonable  to 
assume  that  no  coalescence  of  droplets  occurs  and  that  the  droplet 
size  reached  at  the  throat  will  be  maintained  thereafter  to  the 
nozzle  exit.  That  final  droplet  size  will  be  important  for  the 
performance  of  a two-phase  turbine. 

As  far  as  the  workings  of  the  computer  program  for  the  slip 
variation  are  concerned,  no  peculiarities  or  difficulties  were 
encountered  in  the  expansion  of  a two-phase  mixture  in  the  con- 
verging straight  conical  part  of  the  nozzle.  Much  time,  however, 


was  absorbed  in  interpreting  the  results  for  the  diverging  part 
of  a nozzle.  For  small  nozzles  (with  one  centimeter  inlet  diam- 
eter) the  slip  ratio  was  large  and  the  further  expansion  seemed 
to  require  a departure  from  the  straight  cone  geometry.  A trial 
and  error  adjustment  of  the  cone  angle  was  used  to  arrive  at 
favorable  conditions.  See  the  end  of  Part  I for  further  dis- 
cussions . 

On  the  whole  it  is  felt  that  basic  nozzle  design  tools  were 
developed  in  the  form  of  equations  [for  example  Equations  (21), 

(23),  (26),  (27),  (30),  (34),  (37),  (38),  (57),  (58),  (79),  (82) 

(83),  (86)  and  (90)],  and  in  the  form  of  three  computer  programs: 
for  the  initial  condition,  for  the  limiting  throat  conditions 
and  for  the  study  of  the  entire  expansion.  An  optimum  nozzle 
design  could  not  be  concluded  at  this  point,  since  no  specific 
conditions  or  applications  were  defined  as  a requirement.  More 
work  is  needed  to  define  certain  types  of  nozzle  conditions  and 
the  corresponding  design  recommendations  and  performance  poten- 
tials . 

Since  the  nozzle  performance  study  of  Part  I forms  the  basis  for 
further  two-phase  turbine  work  (Part  II),  the  latter  could  only 
be  advanced  to  the  development  of  basic  procedures  for  its  design 
work.  Once  the  preatomization  and  nozzle  expansion  problem  is 
sufficiently  organized,  more  specific  conclusions  may  be  drawn 
for  the  turbine  performance  potentials  in  terms  of  the  steam 
conditions,  sizes  and  ratings.  Since  development  work  is  in 
progress  with  liquid  turbines  a separate  similar  program  for  two- 
phase  turbines  should  be  initiated. 

For  a further  discussion  of  results  see  "Conclusions  and  Recom-  j 

mendations" . 


PART  I 

ONE  DIMENSIONAL  MODEL  FOR  NOZZLE  WITH  SLIP 


In  the  IHE  model  it  was  assumed  that  the  velocity  of  the  droplets, 

, equals  exactly  the  velocity  of  the  gaseous  carrier,  C&  . 

If  we  allow  now  a slip  ratio 

K , Ca/Cb  > 1.0 

in  a one-dimensional  flow,  (that  is  in  a flow  where  uniform 

* 

conditions  prevail  across  any  section  at  the  length  coordinate  x ) , 
we  must  also  distinguish  between  a volume  or  void  fraction  a of 
the  gaseous  phase  (out  of  the  total  volume  present  at  static  con- 
ditions), and  the  vapor  flow  fraction,  or  quality,  x which  is  the 
ratio  of  the  gaseous  mass  flow  sate,  referred  to  the  total  mass 
flow  per  unit  time  at  a given  section. 

While  the  mass  density  of  the  gaseous  phase,  referred  to  its 
volume  is  1/v  , the  mass  density,  referred  to  the  total  volume 

is  a/v  . The  gaseous  mass  flux  density,  referred  to  the  total 

cl  t # 

cross-section  is  m /A  = aC  /v  . By  definition  m /A  * xmT/ A . 

cL  d cl  cl  1 

Consequently,  the  total  mass-flux  density  is 


ii»T/A 


aC , 


vax 


Cl) 


Considering  that  the  static  volume  fraction  of  the  liquid  is 
1-a  and  its  mass  flow  fraction  1-x,  the  corresponding  equation 
for  the  total  mass-flux  density  is  in  terms  of  v^  is 


(l-ct)C. 

VA  ’ ^Tl^y 


(2) 


9 
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Equating  (1)  and  (2)  allows  a to  be  calculated  when  x , k 
and  the  density  ratio  vb/vfl  ^nown: 


i 


a y = x _JL  rj") 

(l-a)K  (1-x)  vb  1 J 

The  static  void  fraction  a corresponds  only  to  the  density  cor- 
rected flow  fraction  x if  the  slip  ratio  K is  unity;  otherwise 
differences  in  velocity  bring  about  changes  in  flow  rates  for  the 
same  void  fractions.  Equation  (3)  solved  for  1/a  and  l/(l-a)  gives 


1 = 1 ♦ 


:i-*)  Ik  K 


= i + 


TT^T  VT  K 


The  basic  equations  will  now  be  formulated  for  the  momentum, 
energy  and  continuity  of  the  nozzle  flow,  following  Reference  5 
Neglecting  frictional  shear  forces  at  the  nozzle  wall  the 
momentum  equation  for  steady  flow  is 


Adp  = -d(maCa  * mbCb) 


-d[mT(xCa+(l-x)Cb)] 


Since  m^.  is  constant  along  the  nozzle 


-_L  d[xc  ♦ U-xKb] 


The  velocities  may  be  eliminated  using  Equations  (1)  and  (2) : 


Substitution  of  Equations  C4)  and  (5)  for  a and  1-a  gives 


(T)d{(T)[x!va  * xC1'x)[t  * VbK]  * tl-x)!vb]}  (9) 


Following  Reference  5,  the  expression  in  brackets  can  be 
expressed  as  the  product  of  the  quantities  x and  v defined 
as  follows: 


X « xv/T  + Cl-Xl 


XV 

v * jf  * (i-x)vb  /ic  . 


Using  x and  v Equation  (9)  becomes 


-Cl)i j!l  a 

\ A / A 


\ 0 


The  Energy  Equation 


With  i defined  as  stagnation  enthalpy  of  the  mixture,  we  get 
o 


C 2 

V 4a*  -r-41-*) 


The  mixture  enthalpy  is 


i = xi  + (1-x)  i, 
m a b 


Therefore 


i ♦ i xC  * ♦ 
m 2 L a 


d-x)rb>]  . 


The  expression  in  the  brackets  may  be  defined  as  Cg2  • (See 
Equation  (21)  below.)  Again  eliminating  C&  and  by  means  of 
Equations  (1)  and  (2)  we  get 

di  = - [x34  + (l-x)3_^_]j  (16; 

m 2 (\  V a (1  -a) 2 J j 

Using  again  Equations  (4)  and  (5)  we  get 

dim  ’ Td|(s)2[x’va!*  x!fl'x)(!VbI  * ) 


+ x(l-x) 


2 /“VaVb 
\ K 


* V"’)  * U-x)‘vb!]|  . 


Using  the  definition  (11)  for  v and  the  additional  abbreviation 


x = xK  + 


-M. 


r 

(r 


L I 


it  is  seen  that  the  expression  in  the  square  brackets  of  Equa- 
tion (17)  can  be  replaced  by  xv2  . The  energy  equation 


finally  reads 


di  = 

-Id  j 

) * v2| 

m 

2 ( 

\ A 

/ 1 

Integrated  between  the  stagnation  point,  where  mT/A  = 0 and 
a general  point  we  get 


X ■ 


!(?)' 


We  are  now  in  a position  to  formulate  the  deviations  in  velocity 

from  the  ideal  slip  and  loss-free  condition.  When  K=1  , x=l  , 

v=v  so  that  (20)  yields 
m J 


i . i . . fa* 

mo  m 2\  A j 2 


which  is  the  kinetic  energy  of  the  uniformly  accelerated  mixture. 
Note  that  Cg  is  not  necessarily  a result  of  an  ideal  isentropic 
expansion . 

Combining  Equations  (1)  and  (4),  the  mass-flow  density  in  terms 
of  the  gaseous  velocity  is 


C C 

a a 


[xv  + (l-x)v.  K]  v /ic 


Equation  (20)  solved  for  m^/A,  and  with  the  definition  (21)  for 
CE  substituted  gives 


h.  . (_i-V 

A \Ji-J 


This  equation  expresses  the  mass-flux  density  in  terms  of  the  slip 
ratio  K,  the  vapor  mass  fraction  x,  the  specific  volumes  and  the 
enthalpy  change  CE2/2. 

Without  slip  the  predicted  mass-flow  rate  is 


/^t\  = fs 

A IHE  V,n 


The  flow  correction  factor  for  the  influence  of  slip  therefore 
becomes 


C*T>K  . / Vm  \7CE\* 

/ \cs/ 


“T^IHE 


The  above  factor  may  be  applied  to  the  formulas  developed  in  a 
previous  analysis  for  slip-free  conditions.  For  a convergent- 
divergent  nozzle  the  conditions  at  the  throat  are  flow  controlling; 
therefore  the  slip  ratio  K and  the  vapor  mass  fraction  x have  to  be 
known  at  the  throat. 


Special  values  for  the  flow  correction  are  as  follows: 


vr  - J\  ; « yr 

K VVX 

n/x~  = v/iT  ; -^=-  s 1.0 

vVx 

The  result  is  that  at  low  x values  the  presence  of  si 
the  flow  over  the  no-slip  prediction. 

Equating  Equations  (22)  and  (23)  yields 


When  x = 

When  x ■ 


ip 


Considering  that  K = Ca/C^  we  also  get 


Note:  Check  of  equations  for  special  cases: 


1)  x - 0 : 


When  x - 0 there  is  no  vapor;  therefore  intuitively: 


2)  x - 1 : 


Cfl  * 0 and  K = 0 , x -►  00  • 

x - K ; Ca  = Cg 
Cb  * Ce/K 


When  x * 1 there  is  no  liquid;  therefore  C^  * 0,  K -►  «. 
tions  are  found  to  work  even  under  extreme  conditions. 


raises 


(26) 


(27) 


The  equa 
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The  significance  of  Equations  (26)  and  (27)  is  that  the  phase 
velocities  C&  and  are  now  expressed  in  terms  of  slip  ratio  K, 
vapor  fraction  x,  and  the  energy  velocity  C£  that  represents  the 
enthalpy  change  of  the  mixture. 

The  thrust  and  nozzle  efficiency  may  now  be  expressed  similarly. 
Thrust 

Neglecting  off-design  conditions  with  possible  deviations  of  the 
static  pressure  inside  the  nozzle  from  the  back  pressure,  the 
axial  thrust  F is  expressed  as 


F 


iiT[xca  ♦ CX-*)Cb] 

A Z 


(28) 


where  the  velocities  and  and  the  vapor  mass  fraction  x are 
to  be  taken  at  the  nozzle  exit  ("N2") . 

After  substitution  of  Equations  (26)  and  (27) 


F 


(l-x)-L 

n/kJy 


N2 


(29) 


The  expression  in  brackets  is  recognized  as  x defined  in  Equation 
(10);  therefore 


(30) 


Again  the  result  is  expressed  in  terms  of  the  slip  ratio  K at  the 
nozzle  exit  and  the  vapor  mass  fraction  x and  the  energy  velocity 


O 


* 


C£  at  that  station.  The  kinetic  energy  Cg2/2  represents  the  en- 
thalpy change  of  the  mixture,  which  may  differ  from  an  isentropic 
change  on  account  of  friction  losses  inside  the  flow.  The  thrust 
correction  factor  on  account  of  slip  is  (when  m^  is  considered  a 
measured  value) 


(31) 


where  K and  x are  taken  at  the  nozzle  exit  station.  If  the  flow 
correction  factor  (27)  is  included,  the  thrust  correction  factor 
is 


(32) 


Thrust  Coefficient 


We  introduce  the  velocity  C = y/l(i  -i  ) 

J s v o s'm 

an  isentropic  expansion  from  the  stagnation 
thrust  coefficient  may  then  be  defined  as 


which  follows 
condition  i 

mo 


from 

The 


Using  Equation  (30)  we  get 


(33) 
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Efficiency 

An  efficiency  can  be  defined  as 


\CJ2  + W2 


'N 


W2 


xCa2  ♦ (l-x)Cb2 


L 


C 2 


(35) 


Substituting  Equations  (26)  and  (27)  for  C&  and  we  get 

nN  X I * (1'x)  4]  (cf)  (36) 

L 

A 

Considering  definition  (18)  for  x reduces  the  expression  in  the 
brackets  to  unity  (as  it  should  be  for  the  fulfillment  of  the 
energy  Equation  (15)).  Therefore 

O 

ft)  (37) 

The  slip  losses  show  up  only  indirectly  in  the  decrease  of  Cg 
below  Cg  due  to  friction  losses  on  account  of  relative  motion 
f;  between  the  phases. 


A "thrust  efficiency"  amounts  to  taking  the  square  of  the  thrust 
coefficient,  Equation  (33): 


In  Reference  7 relative  slips  are  defined.  Correspondingly,  the 
difference  between  and  may  be  referred  to  the  energy  velocity 
Cg  or  to  a momentum  velocity  C^,  where 


C 2 


C 2 C 2 


(39) 


*T  ' (°E  vl  ) 


N2 


(40) 


Accordingly,  from  Equations  (26)  and  (27) 


V! 


>/kI"  y/i 


r ' yir ) 


(41) 


M 


Ca  - Cb 
CM 


vT  ■ m 


X 


j ' yr ) 


(4  2) 


The  relative  slip  S depends  therefore  not  only  on  K but  also  on 
the  vapor  mass  fraction  x. 

Special  cases: 


0;  x 

/v 

l;  x 


l 

K 


K , x 


1 

s/K 

Vk 


K-l  = S, 


M 


1 - — = S 
1 K 


(43) 

(44) 


19 


The  Determination  of  the  Slip  Ratio  K 


The  above  relationships  allow  the  calculation  of  mass-flow,  thrust 
and  efficiency  when  the  slip  ratio  K,  the  vapor  mass  fraction  x, 
and  the  specific  volumes  v&  and  v^  are  known  at  the  throat  and 
exit  sections. 

In  a first  part  it  will  be  shown  that  the  combination  of  the  momen- 
tum Equation  (12)  with  the  energy  Equation  (19)  yields  an  ordinary 
differential  equation  of  first  order  in  K,  provided  the  losses  are 
known.  This  first  part,  which  excludes  the  loss  formulation,  is 
independent  of  the  form  of  the  flow,  that  is,  whether  one  or  each 
of  the  phases  form  a connected  space. 

To  this  point  the  procedure  will  run  parallel  to  that  of  Reference  5. 
A radical  departure  is  then  taken  in  the  loss  formulation  from  the 
impact  model  of  Reference  5 to  the  frictional  drag  of  a swarm  of 
droplets,  a procedure  more  similar  to  Reference  7. 


I 


Similarly  for  Equation  (19) 


( / m„  \2  . ) /L  V „ . /iHt\  /™t'\ 

|(t)xx  j • * 2 x ''AT/  d(“A~) 


- 2di  = ( 

M*-  -2 

— x v 2 ( 

d (x  v2) 

d(mT/A)  | 

m \ 

A / 1 

x v2 

iftj/A  j 

Substitution  of  Equation  (13)  for  the  mass-flux  density  and  elimi- 
nation of  d (m^/A)  between  Equations  (45)  and  (46)  gives 


m.p/A 


- CF2  d (x  v) 

dim  = * v dp  + -E-  2 

m x 2 / v v 


d (x  v2) 

A ”5 

v yr  *• 


According  to  the  laws  of  differentiation  of  products 


d(x  v) 
x v 


dx  dv 

X V 


d(x  v2)  dx  dv 

A » A A — ’ 

V tf  ‘ V If 


Substitution  of  (48)  and  (49)  into  (47)  yields 


C 2 

* v dp  + -L-  2 

x 2 


dx  dx 

— A 

X X 


_ — 


Substitution  of  Equations  (51)  and  (52)  into  the  last  two  terms 
of  (54)  after  rearranging  gives 


I 


The  first  term  of  Equation  (54)  can  be  rearranged  as  follows 


2x(x  v - nN  x vm) 


dp 

C 2 
LE 


2x(x  v - x vm)  + 2 xx  vm(l  - nw) 


m 


NJ 


dp 

C 2 
LE 


-2xJd -x)(yr  Kvb) 


+ 2 xx  vm(l  - n 


Jp  j 


dp 

C 2 
lE 


(56) 


Recombining  Equations  (55)  and  (56)  yields 


dK 

(x-1) 

dp 

A Y + 

1 '/Va  . rv  \ • V ) 

= 0 

K 

2x (1 -x) 

CLX  ^ 

c 2 
le 

( \K  b/  x(l-x)pK-  -i— j | 

(57) 


Multiplying  by  K/dp  gives  the  alternate  form 


dK  (x-l)K  dx  K ( -xxvm(l  - nM) 


dp  2x(l-x)  dp 


CE2  ( x(l 


:VmU  ~ V + a . Ky 

WB)  (k  b 


0 (58) 


Equation  (57)  or  (58)  is  the  desired  differential  equation  for  K 
with  p as  the  independent  variable.  The  kinetic  energy  C|/2 
represents  the  enthalpy  change  from  the  stagnation  condition  to 
the  pressure  p at  the  inlet  to  the  flow  element  investigated 
between  two  neighboring  sections  inside  the  nozzle.  For  Case  I, 
with  constant  vapor  fraction  x,  the  second  term  for  dx/dp  becomes 
zero.  In  Case  II  with  evaporation,  for  example,  of  wet  steam 
dx/dp  may  be  calculated  according  to  Appendix  II,  which  gives 
relationships  for  the  calculation  of  steam  properties. 
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Solutions  of  Equation  (57)  have  been  obtained  for  an  assumed 

elemental  relative  loss  (1-n^)  for  wet  steam.  It  was  found  that 

certain  limits  in  (1-n  ) have  to  be  observed. 

N 

Solution  of  Equation  (57)  turned  out  to  be  easier  if  (1-n  ) is 

N 

eliminated  by  a certain  loss  model  to  be  treated  next. 

Loss  Model 

Consider  a control  volume  of  length  dx*  and  cross  sectional  area 
A containing  n spherical  droplets  of  diameter  d.  The  ratio  of 
liquid  volume  to  size  of  control  volume  Vc  is 


1 

4> 


Adx* 


(59) 


The  liquid  mass  flow  can  then  be  expressed  with  the  liquid  velocity 
and  its  specific  volume  v^  as  follows 


C.  A 

• _ b 

b ' vb* 


Cbnfe? 

vbdx* 


(60) 


The  total  mass-flow  follows 


mt  = mb/(l“X)  , (61) 

Assuming  a drag  force  Dp  between  droplets  and  vapor  defined  by 
the  drag  coefficient  CQ  and  the  relative  velocity 

D,.  = Cn  1 ■ (Ca'Cb)  tt  d2  . (62) 

F D 

va  2 4 


The  power  loss  of  n droplets  is , neglecting  buoyancy 

L - n Dp(Ca-Cb)  . 


(63) 


. 


-A  : 


» 
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Here  the  droplets  are  considered  spaced  sufficiently  far  apart 
so  that  the  flow  around  the  spheres  is  not  mutually  influenced 
by  neighboring  spheres. 


The  loss  per  unit  mass  follows  by  dividing  Equation  (63)  by 
Equation  (61)  with  mb  according  to  Equation  (60)) 


Tds  = 


nDF (Ca"Cb^  vbdx*(1'x) 


n it  d3  C, 

n F b 


(64) 


Substitution  of  Dp  according  to  (62)  gives 


Tds  - 1Cd(1-xA  d^Cb2(K-l)‘  . 

2 


(65) 


Eliminating  Cb  according  to  Equation  (27),  Cb: 
yields 


C£2/Kx  , 


Tds  = 3 Cn  ^ h?-  (K-  1)’  -L-  . 


a d 


2K 


(66) 


Considering  that  the  relative  loss  is 


1-n 


-Tds 


N 


v dp 
m 


(67) 


we  get 


(68) 
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Accordingly,  the  relative  loss  increases  strongly  with  K ; the 
inverse  dependency  upon  droplet  size  d is  at  first  surprising. 

The  contrary  intuitive  perception  is  that  nozzle  performance 
improves  with  smaller  droplet  size.  It  will  be  seen  that  indeed 
the  final  results  confirm  that  expectation,  after  all  terms  are 
completely  developed,  especially  the  dependence  of  the  slip  ratio 
K on  d.  It  was  later  found  that  Equation  (68)  is  in  agreement 
with  Equation  (78)  and  (79)  of  Reference  9. 

Of  all  the  quantities  occurring  in  Equation  (68)  , the  length  in- 
crement dx*  needs  the  most  elaboration,  since  it  needs  to  be 
expressed  in  terms  of  x , vfl  , v^  , K , , and  dp 
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Calculation  of  the  Diameter  Change  dD  and  the  Length  Increment  dx  . 

One  possible  way  to  express  dx*  is  via  a mass-flux  density  change 
and  therefore  a diameter  change.  If  the  slope  of  the  nozzle  walls 
is  assumed  as  known  by  the  half  angle  6 , the  diameter  change 

dD  is  linked  with  dx*  according  to 

dD  = 2dx*  • tan  6 • (69) 

The  sign  convention  for  the  half-angle  6 can  be  selected  such 
that  for  the  converging  part  of  the  nozzle  6 is  taken  as  negative. 

That  way  dx*  will  be  positive  throughout  the  nozzle. 

The  change  in  mass-flux  density  can  be  linked  to  the  desired 
independent  variables  by  means  of  the  momentum  equation.  Combining 
Equations  (45),  (23)  and  (48)  gives 

-d(m/A)  _ xv  dp  + dx  + dv  (70) 

• — _ _ 

m/A  x Cc2  x v 

E 

{ 

Since  the  mass-flow  rate  is  constant  along  the  axis  the  left-hand 
side  of  Equation  (44)  can  be  reduced  to 


-d (m/A)  = 2dD 

m/A  D 

for  a circular  cross-section  of  diameter  D. 

Using  Equation  (71),  the  momentum  Equation  (70)  becomes 

2dD  m xv  dp  + dx  + dv 

r x cE«  i 

where  dD  may  be  replaced  by  2dx*  tan  6 
Dividing  the  equation  by  dp  gives 


for  the  momentum  equation, 
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(71) 


(72) 


S = 2dx* 

tan5 

xv  1 ^ dx 

+ dv 

2 dD 

1 

(73) 

dp 

D/2 

x CE2  xdp 

vdp 

D dp 

t 

On  the  left  hand  side  dx*/dp  constitutes  the  reciprocal  of  the 
axial  pressure  gradient. 


J 


? 
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In  the  case  of  single  phase  flow  the  change  in  mass-flux  density 
is  easily  obtained  as 

d(m/a)  = d(pC)  = pdC  ♦ Cdp  . (74) 

In  normalized  form,  using  the  specific  volume  v instead  of  the 
density  p we  get 

d (m/A)  dC  dv  (75) 

(m/A)  C v 

Euler's  equation  for  one-dimensional,  gravity-free  isentropic  flow 
yields 

dC  m - vdp  . (76) 

C C2 

Equations  (75)  and  (76)  combined  give 


(77) 


I 

A comparison  of  the  equation  with  Eq.(70)  shows  the  meaning  of 
the  individual  terms  in  Equation  (70)  . xvdp/xC^2  represents  the 
effect  of  a relative  velocity  change;  dx/x  constitutes  the  effect 
of  a change  in  vapor  fraction;  dv/v  represents  the  influence  of 
a specific  volume  change.  For  single  phase  flow,  in  the  converging 
part  of  a nozzle  the  effect  of  the  velocity  increase  more  than 
offsets  the  increase  in  specific  volume,  therefore  the  mass-flux 
density  progressively  increases;  in  the  diverging  part  the  opposite 
is  true.  For  two-phase  flow  the  additional  change  in  vapor  fraction 

i 


•3 


-d(m/A)  _ vdp  + dv 

m/A  Cl  v 


I 


ll 

appears  and  the  effect  of  slip  is  present  in  the  quantities 
t x,  x and  v . To  show  the  effect  of  changes  in  the  slip  ratio, 

dK  , the  differentials  dx  and  dv  have  to  be  expressed  in 
more  basic  terms. 

On  the  right  hand  side  of  Equation  (73)  the  differential  dx  was 
previously  given  in  Equation  (51)  . The  differential  dv  is 
similarly  derived  from  the  definition  of  v in  Equation  (11)  (when 
neglecting  the  change  in  the  specific  volume  of  the  liquid) 

(78) 


xva  r- 1 

dv  * Tyf  '(1'x  Vb^. 

2K 

dx  * W dv’ 

Substitution  of  dx  and  dv  from  Equations  (51)  and  (78)  gives 
for  the  momentum  Equation  (73) 

(79) 


Again,  the  right  hand  side  of  Equation  (79)  is  a function  of  x, 
va » vb,  K,  and  CE 

A discussion  of  the  signs  of  each  term  will  show  whether  its 
influence  is  increasing  the  mass-flux  density  (for  positive  terms) 
or  the  opposite.  In  the  expansion  of  wet  steam  for  example,  the 
last  term  in  Equation  (79),  x dvfl/dp  v/k  is  always  negative;  Z 
is  always  positive,  therefore  Zdx/dp  is  always  negative.  The 
first  term  xv/xC£2  is  always  positive.  Since  the  first  term  in 
V is  negative  for  small  vapor  fractions  x , and  since  the 
second  term  of  V approaches  -1  in  value,  a flow  with  increasing 
slip  (dK/dp<o)  will  tend  to  give  an  increasing  mass-flux  density. 
If,  in  the  overall  effect,  S is  positive,  the  nozzle  will  be  con- 
verging, with  tan6  being  negative. 


The  relative  loss,  Equation  (47)  , may  now  be  reformulated  in  terms 
of  S = 4dx*  tan  6/D  dp  . For  this  purpose  we  define  the  fol- 
lowing : 


W = — C 


(1-x)  vb  (K-l) 


D 


D 


x v. 


d tan6 


(82) 


Equation  (68)  in  terms  of  S and  W now  becomes 


(83) 


Since  W is  dimensionless  and  S has  the  units  of  reciprocal 
pressure,  the  dimensions  check. 


■ -i 


> i 
! 

1 1 
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The  important  conclusion  is  that  the  relative  loss  is  not  only 
determined  by  the  slip  ratio  K and  the  droplet  size  d which  are 
represented  by  W but  the  normalized  inverse  pressure  gradient 
SCf2/^  : good  local  efficiencies  are  possible  with  large  K 
ratios . 
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Droplet  Size 


In  order  to  evaluate  W with  Equation  (82),  the  droplet  size  d 
must  be  known,  so  that  the  drag  coefficient  Cp  can  also  be 
evaluated. 

A criterion  for  the  droplet  size  is  that  of  a limiting  Weber  Number 


We 


(ca  - cb)2d 


<6 


2 v a 
a 


(84) 


where  a is  the  surface  tension,  which  is  given  in  Appendix  II. 
The  droplet  size  limitations  are  given  by 


d < 12  va°  . (85) 

This  criterion  was  used  by  Elliott  in  Reference  7;  it  also  implies 
that  the  droplets  are  spherical  in  shape,  and  that  they  are  all  of 
the  same  diameter.  A further  discussion  of  the  break-up  criteria 
was  already  given  for  example  by  J.O.  Hinze  in  1947  in  Reference  10. 

Using  the  definition  of  the  slip  ratio  K and  Equation  (27)  for 
the  liquid  velocity  in  terms  of  the  energy  velocity  Cg  of 

Equation  (21),  Equation  (85)  reduces  to 


, . 12  vo  Kx 

d a 

CE2 (K-l) 2 


(86) 
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The  drag  coefficient  D is  a function  of  the  Reynolds  number, 
defined  as 


Re  = 


v n 
a a 


Reference  7 uses  Stonecypher ' s least-square  fit  to  Perry’s  tabu- 
lation as  follows: 


In  CD  = 3.271  - 0.8893  In  Re 

♦ 0.03417  (InRe) * (88) 

+ 0.001443  (InRe)*, 


within  the  limits 

0.1  < Re  _<  2xl04  . 

Again,  the  Reynolds  number  may  be  expressed  as  a function  of  the 
basic  variables 


d CF  (K-l) 

Re  = £ • 

vana 


In  conclusion,  W may  now  be  calculated  in  terms  of  x.  v . v.  . 

’ a ’ b 

CE,  K,  and  the  viscosity  na  (see  Appendix  II)  by  substituting 
Equations  (86),  (88),  (89)  into  Equation  (82). 


) 


The  Final  Differential  liquation  for  The  Slip  Ratio  K 

We  are  now  in  a position  to  eliminate  the  unknown  relative  loss 
1-n^  in  Equation  (58)  by  the  use  of  the  relation  (83),  where  S 
has  to  be  expressed  by  Equation  (79).  Fortunately  the  elimination 
of  1-n^  is  possible  without  changing  the  character  of  the  differ- 
ential Equation  (58)  to  find 


Equation  (90)  also  includes  a term  with  dx/dp  , the  change  of 
the  vapor  fraction  with  dp  as  the  expansion  progresses.  The 
calculation  of  dx/dp  for  wet  steam  is  given  in  the  Appendix  II. 

With  the  conditions  T,  x,  v , CP,  K,  a,  n , dx/dp  and  the  diam- 
eter  D known  at  the  inlet  to  a concical  nozzle  element,  it  is 
possible  to  calculate  the  droplet  size  d with  Equation  (86), 
the  Reynolds  number  with  Equation  (89),  the  drag  coefficient  Cp 
with  Equation  (88),  so  that  W follows  from  Equation  (82).  The 
derivative  dK/dp  follows  then  from  Equation  (90) , and  the  relative 
diameter  change  S = 2dD/D  dp  from  Equation  (79).  A problem  is 
posed  by  the  dependence  on  of  dK/dp  , of  d(~l/C^2)  , of 
Re  (-Cgd)  and  of  Cp  , and  therefore  W = f(CE)  . The  depen- 
dence on  CE  is  also  present  in  S according  to  Equation  (79), 


While  the  enthalpy  change,  CE2/2  , from  the  stagnation  point  to 

the  inlet  of  the  element  that  is  being  investigated  may  be  assumed 
as  known,  the  outlet  value  must  be  calculated  for  ready  use  in  the 
next  increment. 

The  static  enthalpy  change  of  the  mixture,  across  the  element 
considered,  equals 


di  = nN  vm  dp  * (93) 

The  local  efficiency  follows  from  the  loss  formula  (83),  whereas 
dp  follows  for  wet  steam  from  Clapeyron's  equation,  if  dT  is 
known 


d£ 

dT 


s -s 
a 


b 


V ■ Vi 

a b 


T(v. 


(94) 


where  rc 


is  the  heat  of  evaporation,  given  in  the  Appendix  II. 


At  first  sight  the  elimination  of  (1  - n^.)  in  Eq . (58)  for  dK/dp 
seems  to  make  an  iterative  procedure  unnecessary.  However,  since 


; ' 


C 2 = 
lE 


-2[(i  - i0)  ♦ 


nNVmdP  J 


(95) 


and  Cg  is  present  in  the  formula  for  droplet  size,  for  Re  and 
therefore  IV,  besides  its  direct  presence  in  Eq.(90)for  dK/dp, 
evaluation  of  dK/dp  at  successive  points  ir.  the  course  of  expan- 
sion over  a nozzle  element  is  affected  by  changes  in  dp  within 
the  element. 


It  is  not  possible  to  eliminate  and  at  the  same  time  preserve 
the  form  of  the  differential  Eq.(90)  that  is  explicitly  solved  for 
dK/dp.  The  following  iterative  procedure  was  therefore  followed. 


The  value  of  C^2  at  the  inlet  to  the  nozzle  element  considered  is 
assumed  as  known,  either  by  being  initially  prescribed  or  by  being 
carried  over  from  the  outlet  of  the  previous  element.  As  outlined 
before  the  values  of  d,  Re,  and  IV  may  then  be  determined. 


In  order  to  eliminate  the  effect  of  C^2  on  IV,  a new  quantity  IV"  is 
defined 


W = -CD(l-x)(K-l)5  vbD 

4vmCE2  32  £2K2va2vmtan5  o 


(96) 


which  does  not  contain  droplet  size  and  C^2  directly.  Similarly, 
the  derivative  S (Eq.  (73))  can  be  expressed  as 


S 


V dK  \ 
TK  / 


+ 


» 
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(9  7) 


The  expression  (1  - n^/C^  in  Eq.(58)  for  dK/dp  can  now  be 
replaced  by 


c 2 

*"E 


W'SCE2 


dK 

3P 


0- 


After  substitution  into  Eq.(58)  we  get 


(98) 


dK 


(x-1)  K dx 
2x  fl -x ) 3p 


+ A MV" 


RCc  2 

k t 


1 


A MV"  WW 


X_ 

2K 


(99) 


where 


A’ 


x x v K 
m 


x(l-x) 


F-yH 


(100) 


The  idea  is  to  use  the  value  of  C 2 at  the  element  entrance  for  a 

E 

first  approximation  of  dK/dp,  according  to  Eq.(99),  where  C£2  occurs 
three  times. 


An  improved  value  of  Cg2  is  obtained  as  follows.  From  Eq . (98) 

1 - „N  ■ w-(r  * v,  «)cE*  * rfcE‘  . 

Multiplying  the  above  equation  by  vmdp  (where  dp  is  the  pressure 
change  from  the  inlet  of  the  element  to  the  particular  station  in- 
side the  element  that  is  being  considered)  and  eliminating  Hv  by 


where  Ai  = - iQ,  the  enthalpy  change  from  the  stagnation  condi- 

tions to  the  inlet  of  the  element  considered. 

It  was  found  that  the  minus  sign  in  front  of  the  square  root  of  the 
discriminant  gives  the  desired  results  for  C^2.  The  result  of  Eq.(99) 
is  used  for  dK/dp  in  Eq.(lOl).  The  new  value  of  C^2  is  then  rein- 
troduced into  Eq.(99),  where  IV"  is  previously  adjusted  for  the 
change  in  C^.  The  recycling  is  stopped  when  initial  and  final  values 
of  Cg2  agree  within  a certain  relative  error  (10“*). 


•1 


I 


B 
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The  differential  Eq.(99)  is  an  ordinary  equation  of  the  first  order. 
Basically,  at  each  temperature  the  properties  v , vb  are  known,  and 
if  the  change  in  the  vapor  fraction  dx/dT  is  predictable,  as  it  is 
(see  Appendix  II),  x is  known  and  dv  /dp,  the  change  of  specific 
volume  of  the  vapor  with  respect  to  the  pressure.  If  the  derivative 
dK/dp  were  only  a function  of  the  pressure,  a simple  integration 
would  do,  but  since  the  derivative  dK/dp  depends  on  the  dependent 

V 

variable,  the  slip  ratio,  K,  a more  elaborate  procedure  is  required. 
(That  dependence  on  K is  present  in  the  factor  K itself,  the  compound 
quantities  x,  x,  v,  Q,  V,  2 and  R,  which  are  all  functions  of  t,  x 
and  K) . 

The  Bulirsch-Stoer  integration  algorithm  was  selected,  which  is  des- 
cribed in  Ref.  Il£jl2.  The  algorithm  was  extensively  tested  with  vary- 
ing step  sizes  in  the  nozzle  program  as  well  as  in  the  solution  of 
a test  case  in  form  of  the  differential  equation  y'  = -y  , which 
has  as  an  exact  solution  the  function  y = Ce  x.  Assuming  as  given 
the  function  value  y = 1.0  (C=l)  at  x = 0,  together  with  the  slope 
y*  = -1.0,  the  Bulirsch-Stoer  algorithm  yielded  the  function  value 
y = 0.006738  for  the  chosen  step  size  of  five  (x  = 5).  That 

• e 

agrees  exactly  within  the  six  displayed  digits  with  y = e '. 

The  algorithm  consists  of  two  parts,  a discretization  part  and  an 
extrapolation  part.  As  the  independent  variable,  (for  example, 
temperature),  is  changed  by  a selected  step  size,  approximate  func- 
tion values  are  determined  in  the  first  discretization  part  accord- 
ing to  a procedure  to  be  described  next.  Starting  values  are  the 
function  value  and  the  derivative  at  the  beginning  of  the  step. 
Approximate  function  values  are  found  at  values  of  the  independent 
variable  that  are  found  by  subdividing  into  equal  increments  the 
total  step,  for  example,  by  values  N[H]  = 1,  2,  3,  4,  6,  8,  12, 

16,  24,  32  and  48,  where  H = 0,  1,  2,  ...,  10.  Figure  1 


% 


» 
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illustrates  the  procedure.  Figure  la  shows  the  way  a first 
approximation  for  the  end  point  of  the  step  is  found.  The  same 
method  is  used  for  obtaining  an  approximate  function  value  at 
the  point  nearest  to  the  starting  point  of  the  step,  as  shown 
in  Figure  lb.  The  same  figure  illustrates  the  procedure  used 
to  move  on  to  the  following  point,  C.  The  derivative  or  tangent  is 
found  at  the  preceding  point,  B,  and  a line  is  drawn  parallel  to  the 
tangent,  again  through  the  preceding  point  A.  The  approximation 
to  the  function  value  at  the  end  of  the  step  is  stored  for  each 
H = 0,  1,  2,  ...»  48  into  a matrix  T[0,H]  . 

Second  part:  Extrapolation  to  zero  step  size.  A rational  function 
is  used  for  that  purpose.  A two-dimensional  matrix  T[J,I]  is 
constructed,  with  the  integer  I varying  from  zero  to  ten  corres- 
ponding to  the  change  in  H from  zero  to  ten.  The  columns  are 
designated  by  the  integers  J.  The  approximate  function  values  at 
the  end  of  the  step,  T[0,H]  form  the  column  designated  as  J =0. 

To  the  left  of  it  (,.T  = -1)  a column  of  zeros  is  placed.  The  suc- 
ceeding columns  for  J = 1,  2,  ...»  9 are  found  recursively  from 
the  formula  (Ref.  13). 


T[  J , I ] 


T[J-1 ,1  + 1 ] 


T[.J  - 1 , 1 + 1 ] - TC-T -1,1] 


F h[I] 

1‘ 

T[J -1,1  + 1] 

- T[J - 1 , I]  ) , 

[hl^I+J  J 

TlJ-1,1+1] 

- TLJ-S,I+lj  | 

(102) 


Such  an  extrapolation  matrix  is  to  be  calculated  for  both  the  final 
value  of  the  slip  ratio  K at  the  end  of  the  step,  and  the  final  value 
of  Cg  2,  which  is  a measure  of  the  enthalpy  change  to  the  end  of  the 
step.  Typical  values  for  the  matrices  are  given  in  Table  I.  The 
last  values,  corresponding  to  T[.T=9,I  = 1]  and  E[J  = 9,I  = 1]  are  finally 


I 


I 

Note:  Two  special  conditions  have  arisen  when  the  convergence 
is  very  rapid:  (a)  The  denominator  T[.I-1,  I + l]  - T [ J-  2 , 1 + 2] 
becomes  zero;  in  that  case  the  extrapolation  may  be  terminated 
at  that  point;  (b)  The  entire  denominator  [ ] J (-1  may  become 
zero;  in  that  case  the  very  large  default  value  for  1/0  of  the 
computer  was  used  for  the  quotient. 


t 


t 


t 


I 
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Table  1. 


Solution 


T[j.l]  MATRIX 


sup  size  o . ?5#c 


I\J 

0 

1 

2 

3 

4 

0 

2 . 693794 

2.693661 

2.693690 

2.693685 

2.693683 

1 

2.493694 

2.693690 

2.693676 

2.693683 

2.693682 

2 

2.693692 

2.693686 

2.693683 

2.693681 

2.693681 

3 

2.693689 

2.693684 

2.693681 

2.693681 

2.693679 

4 

2.693686 

2.693682 

2.693681 

2.693680 

2.693680 

5 

2.693686 

2.693681 

2.693680 

2.693680 

2.693676 

6 

2.693683 

2.693681 

2.693680 

2.693679 

2.693679 

7 

2.693682 

2.693680 

2.693679 

2.693679 

0.000000 

8 

2.693681 

2.693680 

2.693679 

0.000000 

0.000000 

9 

2.693680 

2.693680 

0.000000 

0.000000 

0.000000 

10 

2.693680 

0.000000 

0.000000 

0.000000 

0.000000 

I\J 

5 

6 

7 

8 

9 

0 

2.693680 

2.693681 

2.693680 

2.693680 

2.693678 

1 

2.693681 

2.693680 

2.693680 

2.693680 

li.693t>'9l-K. 

2 

2.693681 

2.693680 

2.693679 

2.693679 

0.000000 

3 

2.693680 

2.693680 

2.693679 

0.000000 

0.000000 

4 

2.693680 

2.693679 

0.000000 

o.nnoouo 

o.oooooo 

5 

2.693679 

0.000000 

0.000000 

0.000000 

o.oooooo 

6 

0.000000 

0.000000 

0.000000 

0.000000 

o.oooooo 

7 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

8 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

9 

0.000000 

0.000000 

0.000000 

0.000000 

o.oooooo 

10 

0.000000 

0.000000 

0.000000 

E [J.l]  MATRIX 

o.oOoooo 

O.UOOOOO 

!\J 

0 

1 

2 

3 

4 

0 

2766.481628 

2766.081376 

2765.991353 

2765.951776 

2765.923368 

1 

2766.181429 

2766 . 008b28 

2765.958896 

2765.926913 

2765.911481 

2 

2766.085426 

2765 .976S26 

2765.935236 

2765.913948 

276S. 917195 

3 

2766.037781 

2T05. 950025 

2765.920021 

27oS. 917097 

276S. 922832 

4 

2765.989027 

2765.931152 

2765.917496 

2765.970505 

2705.898256 

5 

2765.963706 

2765.921882 

2765.852380 

2765.900079 

2765.861309 

6 

2765.940470 

2765.894940 

2765.900810 

2765.891655 

2765.874392 

7 

2765.920550 

2765.899610 

2766.008755 

2765.875619 

0.000000 

8 

2765.908917 

2765.857237 

2765 ,877608 

0.000000 

0.000000 

9 

2‘65. 886307 

2765.875760 

0.000000 

0.000000 

0.000000 

10 

2765.880448 

0.000000 

0.000000 

0.000000 

0.000000 

I\J 

S 

6 

7 

8 

9 

0 

2765.910613 

2765.916724 

2765.907135 

2765.896062 

276S. 921942 

1 

2765.916993 

2765.908338 

2765 . 8962SS 

2765.925809 

12765. 8-3004l»CrS 

2 

2765.911894 

2765.896615 

2765 . 939S63 

2765.873445 

0 . oooooo 

3 

2765 . 8971 7S 

2766.027232 

2765.873896 

0.000000 

0.000000 

4 

2765 . 751926 

276S. 874370 

0.000000 

0.000000 

0.000000 

S 

2765.874602 

0.000000 

0.0001)00 

0.000000 

0.000000 

6 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

7 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

8 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

9 

o.oooooo 

o.oooooo 

0.000000 

0.000000 

0.000000 

10 

O.OOOOOO 

o.oooooo 

o.oooooo 

0.000000 

0.000000 

I 
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The  Stepping  Procedure  of  the  Program 

For  wet  steam  as  a medium,  temperature  was  selected  as  the  inde- 
pendent variable.  A step  size  is  selected,  an  initial  nozzle 
diameter,  the  initial  slip  ratio  K is  taken  over  from  the  result 
(final  value)  of  the  previous  step  (concerning  the  value  at  the 
inlet,  see  the  next  section)  together  with  the  vapor  fraction  x 
and  the  enthalpy  change  (C£2/2)  from  the  stagnation  condition 
as  reference.  At  the  beginning  of  each  step  the  mass-flux  den- 
sity m/A  = C / (v  >fx  ) is  calculated.  Since  the  mass-flux  is 
L 

the  same  along  the  nozzle  axis,  the  ratio  of  the  mass-flux  densi- 
ties from  step  to  step  equals  the  inverse  of  the  cross-sectional 
area  ratio,  or 


D 

2 

rr 


I (m/A) t 

\ wrr 


(103) 


A local  efficiency  can  also  be  calculated,  from  the  ratio  of 
enthalpy  change  across  a step  to  vmdp : 


'N 


ACe2/2 


(104) 


The  nozzle  throat  is  defined  by  the  maxiumum  mass-flux  density. 
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The  Initial  Condition 


An  initial  kinetic  energy  can  be  prescribed  together  with  an  inlet 
diameter  of  the  nozzle  and  the  temperature  and  vapor  fraction  con- 
ditions. The  initial  slip  ratio  K can  also  be  selected;  however, 
experience  with  running  the  program  suggests  that  K is  best  se- 
lected such  as  to  make  the  derivative  dK/dp  at  the  inlet  equal 
to  zero,  in  order  to  avoid  sharp  adjustments  in  K at  the  begin- 
ning of  the  expansion  in  the  nozzle. 


When  Eq.  (90)  for  dK/dp  is  set  equal  to  zero,  a condition  for 


a limiting  value  of  is  obtained  as  follows 


where 


W, 


{ep-  (x  * Kvb)  ' ifrxrr  U } Q 


x x R 


+ 

lv 

q = 4xd-x)(yr  - 


R = Z 


dx 


3?  + 


dv 

x a 
Xp" 


}Jk 


yr  - 4= 


yr  + yr 


- vbyiT 


(105) 


(106) 


(107) 


(108) 


The  value  WL  is  a function  of  temperature  (pressure),  vapor 
fraction  x,  v and  vK , dx/dp,  dv  /dp,  K and  the  initially  assumed 

cl  D o 

Cg2.  From  the  definition  Eq.  (82)  for  W follows 


dL  = 3 D vb(l-x)  (K-l)j 

CT  = 2W.tanS  v 51  K 
D La 


(109) 
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In  order  to  solve  for  the  limiting  droplet  size,  dL>  required 
to  make  dK/dp  = 0,  the  expression  for  the  Reynolds  number  (87) 
is  divided  by  Cp  as  follows 


Re  . 


D var|f 


(110) 


Eq.  (109)  for  dL/CD  may  be  substituted  into  the  above  Eq.  (no)  for 
Re/Cp.  Since  the  drag  coefficient  is  a function  of  Re, 


In  CT 


A + By  + Cy2  + Dy 


where 


In  Re  , 


(111) 


(112) 


CD/Re  may  be  expressed  as  a function  of  Re  as  follows 


'(r?)  " A * 


( B - 1 ) y + Cy  2 + Dy  3 , 


(113) 


which  leads  to  a cubic  equation  for  y when  Cp/Re  is  given 


y3  + g y2  + Y y + ,ln(Cp/Re)  = 0 .(114) 


Substitution  of  y - x • leads  to  the  cubic  equation 


where 


x3  + ax  + b » 0 


_ B - 1 1 /C\2 

" ~tt  * t Id) 


(115) 


(116) 


b = A ♦ ln(Cn/Re) 


’e)  2 ( C V 1 C / B-l  \ 

it  U/  ‘ t d \~ir)  * 


The  discriminant  follows 


* • W*  (*)* 


(117) 


If  A < 0 , which  is  usually  the  case,  there  are  three  real 
solutions,  including 

x = 2 cos(4>  + 240°)  . ( 1 ] 8) 

where  <j>  follows  from 


cos  3<J> 


OI9) 


The  values  of  the  constants  are,  according  to  Eq.  (88)  A = 3.271, 
B = -0.8893,  C = 0.03417,  D = 0.001443,  so  that  a = -1496.198  and 
(a/3)’  = -1.240518  x 108  . The  constant  b/2  becomes 


b/2  = 6792.464  + 346.5  ln(Re/CD)  (120) 

and  y = x - 7.89328. 

As  the  slip  ratio  K is  varied,  the  limiting  d^  required  to  make 
dK/dp  = 0 may  be  calculated.  A typical  plot  is  shown  in  Figure  2. 
The  relative  loss  (1  - n^-)  encountered  due  to  the  dragging  of  the 
droplets  may  be  calculated  and  plotted  on  the  same  figure.  At 
the  same  time  the  Weber  Number  criterion  yields  a droplet  size  d^ 
from  Eq . (86),  which  may  also  be  plotted  as  a function  of  K.  It 
is  seen  in  Figure  2 that  the  d^  and  d^,  lines  intersect  at  some  K 
value.  That  value  of  the  slip  ratio  is  the  one  to  be  used  ini- 
tially if  no  pre -atomizat ion  is  provided  before  the  nozzle  and  the 
"Weber  breakup"  is  assumed.  It  was  found  that  indeed  when  such 
an  initial  K value  was  used  with  the  main  program  the  slip  would 
initially  obey  dK/dp  = 0. 


Similarity  - The  initial  value  problem  for  K just  described  gives 
the  possibility  of  studying  the  performance  of  geometrically  similar 
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nozzles  of  a certain  ratio  of  droplet  size  to  nozzle  diameter  at 
the  inlet.  Figures  2,  3,  4 and  5 show  the  "limiting  droplet  size 
d^"  required  to  give  dK/dp  = 0 for  different  size  nozzles  (P  = 
0.005,  0.01,  0.1  and  0.2  m) . While  the  magnitude  of  W^  is 
independent  of  nozzle  diameter,  d^/C^  contains  D in  Eq.(109). 
The  loss  curve  (1  - 0^),  however, is  independent  of  diameter 
since  S of  Eq.  (79)  is  independent  of  D and 
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Similarly,  the  droplet  size  according  to  the  Weber  criterion  is 
independent  of  nozzle  diameter  D,  according  to  Eq . (86) 


12  v oK  x 
j < a 

dW  " ” 2 


ce2(k  - l)2 


= f(t,  x,  K,  C 2) 


(121) 


The  points  of  intersection  of  the  lines  d^  and  d^.  move  to  smaller 
K ratios  as  the  nozzle  size  increases.  The  droplet  size,  follow- 
ing the  dj,r  curve  increases  correspondingly.  The  ratio  D/dE 
varies  with  K and  D as  follows  : 
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The  Throat  Conditions 
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The  throat  is  located  at  the  plane  of  maximum  mass-flux  density 
m/A.  At  the  same  time  the  local  relative  losses  should  reach  their 
minimum  at  the  throat.  Since  the  relative  loss  1 - n^.  is  propor- 
tional to  the  product  W'S  , the  slope  at  the  wall,  tan  6,  cancels 
out,  so  that  a discontinuity  can  be  tolerated  by  the  one -d  imens  ional 
theory.  Since  the  axial  length  increment  dx*  is  theoretically 
zero  at  a sharp  throat,  the  loss  should  accordingly  be  zero.  That 
requires  the  pressure  gradient  dp/dx*  to  tend  toward  infinity, 
see  Eq . (73) . 

In  an  actual  program  execution  it  was  initially  found  difficult  to 
maintain  the  accuracy  required  to  achieve  coincidence  atone  and  the 
same  temperature  of  the  two  events  of  maximum  mass-flux  density  on  one 
hand,  maximum  gradient  and  minimum  efficiency  on  the  other  hand. 

The  difficulty  was  attributed  to  the  effect  of  CE2/2,  (that  is, 
the  value  by  which  the  local  enthalpy  differs  from  the  value  at 
stagnation  conditions).  The  effect  of  is  explicitly  present 
in  the  differential  Eq.  (90),  as  well  as  implicitly  through  the  pre- 
sence in  W,  see  Eq . (96). 


A refined  equation  for  dK/dp  was  first  developed,  based  upon  the 
assumption  of  a negligible  effect  on  droplet  size  of  a varying 
relative  velocity  at  the  droplet  over  the  axial  extent  of  an  ele- 
ment. Next  the  influence  of  a varying  CP  was  considered  every- 
where else  and  was  algebraically  eliminated.  Further  details  were 
already  given  in  connection  with  the  quadratic  Eq . (101)  for  . 


The  question  is  whether  the  peak  in  mass-flux  density 
accurately  determined  when  m/A  depends  on  W,  and  the 
should  be  adjusted  according  to  its  definition  - Eq . 
changing  tan  6 at  the  throat  (from  the  value  in  the 
part  to  the  value  in  the  diverging  part). 


m/A  can  be 
value  of  IV 
(82) , to  the 
converging  • 
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The  answer  is  that  if  (1  - n^.)  -*■  0 at  the  throat,  Eq . (58)  for 
dK/dp  may  be  resorted  to,  which  is  entirely  general,  (the  loss 
mechanism  being  represented  solely  by  (1  - n^)  in  one  of  the 
terms).  Since  a change  in  tan  5 at  the  throat  only  affects  the 
loss  formulation,  a negligible  loss  situation  is  made  independent 
of  5 and  is  represented  by  Eq . (58)  with  (1  - n^.)  -*-0,  that  is 
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(122) 


For  lower  values  of  x and  moderate  K values  we  have  (l-x)  > 0,  and 
dx/dp  < 0.  The  third  term  is  always  positive  when  va>>vb»  ar>d  much 
larger  than  the  second  term  of  Eq.(122).  As  a result  dK/dp  is  gen- 
erally negative,  that  is,  the  slip  increases  at  the  throat.  That 
is  borne  out  by  actual  computer  runs. 


The  other  relation  that  has  to  be  met  at  the  throat  is  Eq . (73), 
which  represents  the  reciprocal  of  the  normalized  pressure  gradient. 
With  the  pressure  gradient  tending  towards  infinity,  S becomes  zero. 

If  the  temperature  (or  pressure)  and  the  vapor  fraction  is  assumed 
as  known  at  the  throat,  the  enthalpy  difference  C^2/2  may  be  cal- 
culated from  Eq.(73)  for  an  assumed  slip  ratio  K as  follows: 


where 
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were  given  before. 

A check  may  be  obtained  for  Eq . (123)  by  setting  conditions  for 
single  phase  operation  with  an  ideal  gas  as  follows: 

For  K = 1.0  (no  slip)  : 

x-*-1.0  , x + 1.0  , v -*•  v 
’ ’ m 


V = 2x  - 1 - [xva  - (l-x)vb]/vm 


For  an  Ideal  ('as: 


so  that 


dv  _ - 1 dp 

v X p 


After  substitution  into  Eq . (123)  we  get 


C 


E 


2 


kpva 


(124) 


That  is  the  square  of  the  sonic  velocity  in  an  ideal  gas. 


Returning  now  to  two-phase  conditions,  Eq.  (123)  may  be  programmed 
and  results  obtained  as  presented  in  Figure  6 for  the  variation 
of  with  changing  slip  ratio  K,  and  throat  temperature  as  a 
parameter.  Additional  exploration  over  a wide  range  of  K values 
shows  that  for  a given  temperature  reaches  a minimum  value  for 
a certain  slip  ratio.  Over  a range  of  temperatures  from  100  to 
200°C,  that  minimum  value  changes  only  from  about  118  m/s  to  120  m/s. 

Knowledge  of  the  value  of  and  the  slip  ratio  K at  the  throat 
(for  given  temperature  and  vapor  mass  fraction  x)  allows  the  cal- 
culation  of  the  mass-flux  density  m/A  at  the  throat,  which  is 
useful  for  the  sizing  of  the  throat  for  a required  mass  flow. 

Typical  mass-flux  densities  are  shown  in  Figure  7 for  x = 0.125 
and  t = 100  to  150°C  at  the  throat. 
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In  Reference  4 a normalized  mass-flux  density 


0 = m J RTq  /ATpo 


(125) 


was  presented  based  upon  the  slip-free  IHE  model,  and  a given 
nozzle  inlet  temperature  in  the  range  from  100°C  to  300°C. 

Note  that  according  to  Equation  (25)  an  increased  mass-flux 
density  at  the  throat  is  due  to  the  ratio  vm/v  /k7  which  tends 
toward  K for  small  x,  and  the  speed  ratio  (CE/Cg)*. 

It  is  necessary  to  distinguish  clearly  between  three  situations: 

(1)  the  conditions  and  results  of  the  IHE  model,  (2)  the 
results  of  the  slip  model  at  K>1,  (3)  the  results  of  the  slip 
model  at  K=l.  A comparison  of  mass-flux  density  may  be  made 
based  upon  either  (a)  the  same  inlet  stagnation  conditions  or  (b) the 
same  pressure  and  temperature  condition  at  the  throat.  Equation 
(25)  yields  for  K=1  (no  slip) 

(mt)  K=1  = / £e\* 

(lV  IHE  \Cs/ 


If  the  temperature  at  the  throat  and  at  the  inlet  is  assumed  the 
same  for  both  models,  the  enthalpy  drop  should  also  be  the  same 
for  no  slip,  that  is  C^*  = Cg*  and  the  results  of  the  two  models 
(1)  and  (3)  should  coincide  for  K=l. 

Since  Equation  (123)  gives  a value  of  C£2  which  is  only  de- 
pendent on  the  quantities  t,  x and  K,  it  is  not  tied  to  a 
certain  inlet  condition  tQ,  except  that  it  should  not  exceed  Cg 
and  therefore  calls  for  a maximum  tQ.  The  above  assumption  of 
equal  stagnation  conditions  for  the  two  models  when  the  throat 
condition  is  the  same  is  therefore  not  necessarily  justified,  and 
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larger  mass-flux  densities  are  possible  for  the  slip  model  when 
the  throat  conditions  are  reached  at  a lower  enthalpy  than  for 
the  I HE  model  (starting  from  the  same  stagnation  enthalpy) . 

The  curves  for  Cp  of  Figure  6 show  a strong  increase  of  Cp  towards 
lower  slip  ratios  K.  Its  influence  on  mass-flux  density  seems  to 
override  the  effect  of  the  first  factor  of  Equation  (25),  since 
Figure  7 shows  the  highest  mass-flux  density  at  small  slip  ratios. 
Since  Cp  is  only  weekly  dependent  upon  temperature,  according  to 
figure  6,  the  question  of  how  the  mass-flux  densities  compare  for 
the  two  models,  based  upon  identical  stagnation  conditions  may 
have  to  be  answered  by  actual  computer  runs  for  both,  the  IHE 
model  and  the  slip  model  (in  order  to  locate  the  throat). 

From  a series  of  such  computer  runs  it  was  borne  out  that  the 
limiting  values  of  Cp  and  m/A  given  in  Figures  6 and  7 are 
closely  approached  if  the  local  efficiency  is  high.  Low  local 
efficiencies  near  the  throat  are,  however,  possible  if  the  pres- 
sure gradient  is  low.  Throat  values  of  C£  and  m/A  are  then  much 
lower  than  the  limiting  values  given. 

A basic  question  touched  upon  is  that  of  the  value  of  the  sonic 
velocity  at  the  throat  in  a two-phase  mixture  with  slip.  For 
further  details  see  References  6,  9 and  14 

The  limiting  value  of  Cp  given  as  a function  of  K allows  also  the 
determination  of  a minimum  droplet  sire  as  a function  of  K that 
may  be  expected  at  the  throat  and  thereafter.  Figure  8 shows  such 
results.  It  is  seen  that  values  of  d of  the  order  of  a few 
microns  are  possible  if  high  local  efficiencies  near  the  throat 
are  realized. 
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Figure  8.  Throat  Condition:  Minimum  Droplet  Size 


Experience  With  Running  The  Computer  Program  for  the  Nozzle 
Expansion 

The  program  was  systematically  debugged;  the  results  of  the  sub- 
routine for  the  properties  of  wet  steam  were  checked  against 
tabular  values,  until  the  deviations  were  a fraction  of  one  per 
cent.  For  the  viscosity  of  saturated  steam  the  latest  interna- 
tional standards  were  initially  used;  in  order  to  reduce  the 
program  running  time  a formula  of  the  form  of  SUTHERLAND'S  was 
selected  and  the  two  constants  adjusted  such  that  close  agree- 
ment with  the  "International  Formula"  was  achieved,  see 
Appendix  II. 

Most  running  experience  was  obtained  with  150°C  inlet  temperature, 
an  initial  vapor  mass  fraction  of  x=  0 . 1 , an  initial  velocity  CE  of 
17.4  m/s,  and  a very  small  nozzle  inlet  diameter  of  0.01  m (see 
Figure  9).  Without  preatomization,  the  Weber  Number  criterion 
gives  an  initial  droplet  size  (for  dK/dp=0)  of  190  pm,  and  an 
initial  slip  ratio  of  K=4.30=Ca/C^  . It  climbs  initially  when 
using  6=  - 15°  =constant  over  the  converging  part.  The  slip  ratio 
thereafter  decreases  and  reaches  a minimum  of  K=4.05  at  about 
143°C  and  then  increases  again  towards  the  throat,  where  it 
reaches  a value  of  K=5.70  at  a temperature  of  126. 5°C.  The  vapor 
fraction  at  the  throat  is  0.1399,  the  mass-flux  density 
m/A=3202.2  kg/sm2 , the  droplet  size  d=11.6  pm,  the  velocity 
corresponding  to  the  enthalpy  change  Cj:=141.7  m/s,  the  throat 
diameter  D=5.702  mm  and  an  S-value  of  S=1.60  x 10'®  m2/N.  A 
sharp  transition  from  one  cone  angle  to  another  of  opposite  sign 
is  permissible  at  the  throat;  the  K curve  will  then  show  a dis- 
continuity of  the  derivative,  whereas  the  efficiency  and  the  gra- 
dient S will  experience  step  changes.  As  the  angle  6 changes 
its  sign,  W will  also  change  sign  from  positive  to  negative,  and 
S will  change  sign  from  plus  to  minus. 

While  the  normalized  inverse  of  the  pressure  gradient,  S,  may  be 
used  to  calculate  the  diameter  change  for  a certain  temperature 
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or  pressure  change,  it  is  more  accurate  to  calculate  the  mass- 
flux  density  from  Equation  (23)  at  each  inlet  to  an  elementary 
step,  and  use  its  ratio  across  a step  to  determine  the  area  ratio 
required.  Still,  the  quantity  S gives  an  indication  of  the 
diameter  changes  and  the  local  efficiency  potential. 

Near  the  throat,  S tends  to  be  very  small  and  the  pressure 
changes  are  large  for  a given  temperature  change. 

In  the  diverging  part  of  the  nozzle,  a negative  S should  increase 
in  its  absolute  value  in  order  to  allow  the  expansion  to  proceed 
with  smaller  and  smaller  pressure  changes  (for  given  uniform 
temperature  changes).  At  the  same  time  the  losses  are  bound  to 
increase  according  to  Equation  (83) . 

For  very  small  nozzles  (of  the  order  of  0.5  cm  throat  diameter) 
the  slip  ratio  was  large,  and  orderly  expansion  in  the  diverging 
part  was  found  difficult  to  obtain  with  a simple  straight  cone 
geometry.  The  program,  however,  allows  changing  the  cone  angle 
6 at  every  elementary  step  such  that  a compromise  in  the  varia- 
tion of  S and  is  obtained. 

N 

With  a prescribed  fixed  cone  angle  in  the  diverging  part  of  the 
nozzle,  the  situation  was  encountered  where  a negative  S would  in- 
crease in  absolute  value  with  progressing  expansion.  While  the 
mass-flux  density  was  still  decreasing  from  step  to  step,  the 
pressure  changes  were  large,  the  efficiency  was  rising,  and  even- 
tually, when  S was  allowed  to  change  sign  from  the  normal  negative 
to  a positive  sign,  the  pressure  gradient  would  reverse  and  the 
efficiency  would  reach  unity  at  the  same  point.  Within  the  con- 
straints of  the  prescribed  equation  system,  which  does  not  make 
allowance  for  the  occurrence  of  shocks,  only  an  increase  in  cone 
angle  6 downstream  of  that  point  could  avert  the  physically  un- 
tenable condition  of  a rising  pressure  linked  to  a prescribed 
falling  temperature.  The  limiting  condition  as  described. 
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gradually  develops  when  the  absolute  value  of  S has  been  de- 
creasing steadily,  which  means  that  for  uniform  steps  of  tem- 
perature the  diameter  increase  and  therefore  the  axial  movement 
of  the  state  point  along  the  nozzle  has  become  less  and  less. 

That  condition  benefits  the  local  efficiency  because  the  drag 
work  becomes  less  and  less,  until  (when  S changes  sign  and  goes 
through  zero)  the  losses  go  to  zero  when  the  axial  displacement 
dx*  has  ceased  to  develop.  Obviously,  at  this  point  no  further 
regular  expansion  is  possible:  the  selected  cone  angle  6 is  in- 
sufficient to  accommodate  an  increasing  volume  flow.  The  above 
development  goes  parallel  with  a steadily  increasing  slip  ratio  K. 


A guideline  for  the  required  cone  angle  6 may  be  obtained  by 
solving  for  W while  dK/dp  of  Equation  (90)  has  a prescribed  in- 
creasing value  (dK<0) . The  equation  for  becomes  accordingly 
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The  cone  angle  6 required  follows  from  Equation  (82).  Use  of 
Equation  (126)  in  the  program  leads  to  an  automatic  calculation  of 
the  required  nozzle  geometry  that  will  give  a prescribed  slip 
ratio  slope  at  the  inlet  to  each  nozzle  element. 


Such  a sample  solution  is  given  in  Figure  10. 
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PART  II 


DROPLET  TRAJECTORIES  IN  A TURBINE  PASSAGE 

In  a cylindrical  coordinate  system  (r,0,z)  the  material  derivative 
is  (Cr  = radial,  Cy  = peripheral  velocity  components  of  the  droplet) 


_D_  = J_  - r _L  + _Ji  _L  + r * 

Dt  3 1 r3r  r 30  z3z 
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The  projections  ar»au>a-7  of  the  acceleration  vector  of  a liquid 
droplet  in  the  coordinate  directions  are 
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If  we  neglect  the  force  due  to  the  pressure  gradient  in  the  passage, 
the  buoyancy  and  mass  changes  of  the  droplet,  the  aerodynamic  drag 
is  the  only  remaining  force  with  the  drag  coefficient  Cn 


D = CDApa 


CC 


(129) 


Using  the  mass  of  a spherical  droplet 
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and  equating  D to  ma  we  get 
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(131) 


Equating  the  individual  components  we  obtain 
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Such  a set  of  equations  was  solved  for  the  two-dimensional  case 
(C,  = 0)  in  Reference  15  by  numerical  integration.  Unfortunately 
insufficient  informat  ion  is  g i ven  there  for  generalizing  the 
results . 


A simple  analysis  of  droplet  trajectories  was  developed  in  Re- 
ference 4.  The  result  will  be  given  by  the  ratio  of  the  radial 
drift  velocity  w^  of  the  droplets  in  relation  to  the  steam  ve- 
locity C.  A flow  angle  deviation  of  the  droplets  from  the  steam 
is  thereby  given. 

Since  considerable  analysis  was  applied  to  the  study  of  the  processes 
in  the  nozzle  (because  they  form  the  foundation  of  what  is  to  happen 
thereafter  in  the  turbine)  the  criteria  developed  in  Reference  4 
will  be  repeated  here,  together  with  some  sample  calculations. 

The  relative  velocity  between  droplet  and  steam  in  radial  direction 
is  designated  as  wb  and  needs  to  be  calculated. 
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The  drag  force  on  the  droplet  of  diameter  d is  defined  by 


D 


(133) 


where  v is  the  specific  volume  of  the  surrounding  medium.  The 
drag  coefficient  Cp  is  given  as  a function  of  Reynolds  number 
Re  = w,  d/v  in  Figure  11  according  to  a relation  given  in  Re- 
ference  7 which  is  valid  for  0 . 1 < Re  < 20,000:  [See  also  Fq . (88).] 

In  CD  = 3.271  - 0.8893  In  Re  + 0.03417(ln  Re)1 2  + 0.001443(ln  Re)3 

Since  the  slip  velocity  w^  is  not  initially  known  for  the  calcula- 
tion of  Re  and  Cp,  the  quantity  Re^/C^  was  plotted  against  Re  in 
Figure  12  , since  Re^/C^  is  expressible  in  terms  of  known  values 
as  follows. 

The  inertia  force  on  a spherical  droplet  of  diameter  d and  speci- 
fic volume  v^  that  travels  with  velocity  C on  an  initially  assumed 
circular  path  of  radius  of  curvature,  r^,  is 


1 ird 3 C 2 
Vb  ^ 


Corrected  for  buoyance  B we  get 
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Equating  the  two  forces  of  Equations  (133)  and  (135)  gives 
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It  is  convenient  to  refer  velocities  to  the  reference  velocity 


e = J P v and  lengths  to  the  mean  free  path  length 


Accordingly 
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(139) 


Therefore,  the  proof  is  furnished  that  the  quantity  Re>/c^~  is 
solely  composed  of  known  quantities,  namely  the  droplet  diameter 
d,  the  vapor  property  Xa,  the  ratio  of  specific  volumes  v^/v^, 
the  radius  of  curvature  of  the  path,  r^,  and  the  normalised  abso- 
lute velocity  C/e  ('v  Mach  number)  along  the  path.  Entering 
Figure  12  on  the  ordinate  will  yield  Re  and  consequently  the 
desired  w^/e  from  (138) 
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Equation  (139)  shows  that  the  expression  Re>/c^  is  a function  of 
the  relative  droplet  velocity  C/e  in  the  path  with  radius  of 
The  resulting  relative  slip  velocity  w^/e  is  a 


curvature  r 


r-f 


(J 


measure  of  the  deviation  from  the  original  path.  For  this 
purpose  the  ratio  w,/C  is  significant  and  may  be  used  for  a 
step  by  step  trajectory  analysis. 

An  illustrative  example  might  be  helpful.  Assume  an  absolute 
pressure  of  steam  of  0.250  bar  that  corresponds  to  a saturation 
temperature  of  65°C.  Accordingly 


6.206 


10 


-3 


= 6206 


vj p v = 393.9  m/ s 
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vana  _ 6.206  x 1.065x10  5 , c-,a  in-7 

e 39379 1.678  x 10  m 


d = 2 pm 

« 11.919 

a 

r^  = 10  cm 

C = 250  m/s 

| = 0.635 


0.1  m 


Equation  (139)  yields 
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"b  _ 0.0386 


or  an  angle  deviation  of  3.5  . 


0.0608 


If  the  droplet  size  is  diminished  to  1 pm  the  result  is 


1.0885 


= 0.06 


wb/e 


0.06 

37Tb 


= 0.0101 


wb/C  = 0.0159 


a = 0.908 


It  is  seen  that  the  deviation  of  the  liquid  droplets  from  the 
path  of  the  steam  due  to  drift  is  small  for  such  small  droplet 
sizes  and  a radius  of  curvature  of  10  cm,  see  Figure  13. 

To  summarize,  the  non-dimensional  parameters  which  control  the 
droplet  trajectories  in  the  turbine  are 


— the  Mach  index  of  the  steam  flow  , 


d the  ratio  of  droplet  size  to  the  radius  of 
r^  curvature  of  the  path  in  the  bucket  , 


d the  ratio  of  droplet  diameter  to  the  mean 

X free  path  length  of  the  molecules  in  the  steam, 

3 


— the  specific  volume  ratio  of  steam  to  liquid. 
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The  smaller  Re^/C^  , 

the 

smaller  Re,  and  therefore 
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The  lower  the  specific  volume  va  of  the  steam,  the  lower  the 
let  deviation  w^/e.  For  Re  = 0.1  the  simple  Stokes  formula 
valid 

drop  - 

is 
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Therefore , 

Res/C^ 

= s/24  Re 

and  after  substitution  into  Equation  (139)  we  get 
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The  normalized  deviation  w^/C  in  the  Stokes  regime  is  propor- 
tional to  droplet  diameter  squared,  the  Mach  index 


C/  si p v. 


the  ratio  v /v^  (if  va>>vb^  anc*  inversely  proportional  to  the 
product  rjc^a*  For  further  discussion  s'ee  ” Conclus ions  and 
Recommendations” . 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Two  areas  of  two-phase  flow  were  considered:  (a)  the  one- 
dimensional flow  of  wet  steam  in  a converging -diverging  nozzle, 
(b)  the  two-dimensional  flow  through  impulse  turbine  buckets. 

A theory  based  upon  a simple  model  was  developed  for  the  nozzle 
problem,  which  is  a foundation  for  the  turbine  problem. 

Both  the  slip  ratio  K and  the  droplet  size  d at  the  nozzle 
discharge  are  needed  for  the  solution  of  the  turbine  problem. 

To  determine  the  slip  ratio,  a proven  local  relative-loss  model 
was  used  to  express  the  frictional  drag  losses  between  droplets 
and  steam,  while  the  droplet  size  was  established  by  limiting 
the  value  of  the  Weber  number  to  six. 

Expressions  for  the  mass-flux  density,  vapor  and  liquid  velo- 
cities, droplet  size,  thrust  and  overall  nozzle  efficiency  have 
been  included.  These  relationships  have  been  solved  numerically 
by  a digital  computer. 

The  results  predict  an  increased  mass-flux  density  is  to  be 
expected  in  the  presence  of  slip  for  mixtures  of  low  vapor  mass 
fractions  compared  to  the  prediction  of  the  IHE  model.  In  addi- 
tion, the  model  indicates  that  larger  slip  ratios  do  not  neces- 
sarily mean  large  losses  and,  hence,  low  local  efficiencies, 
although  this  is  generally  the  case  for  the  overall  efficiency. 

In  the  absence  of  preatomization,  purely  algebraic  considera- 
tions establish  the  initial  slip  ratio  and  droplet  size  at  the 
nozzle  inlet.  If,  on  the  other  hand,  the  droplet  size  is  pre- 
determined by  means  preceding  the  nozzle,  the  corresponding 
initial  slip  ratio  can  also  be  determined  algebraically. 

The  ratio  of  the  nozzle  and  droplet  diameters  is  predicted  to 
increase  as  the  nozzle  diameter  increases,  which  results  in  a 
higher  nozzle  performance.  For  a given  nozzle  size,  smaller 
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droplets  (below  the  Weber  threshold)  result  in  lower  values 
of  slip  and  local  losses  if  K < 2.5. 


In  the  development  of  the  computer  program  for  the  step-wise 
nozzle  expansion,  entirely  different  situations  were  encountered 
in  the  converging  as  compared  to  the  diverging  part  of  the 
nozzle.  While  the  assumption  of  a given  cone  angle  throughout 
the  converging  portion  lead  to  solutions  without  difficulties, 
a similar  assumption  for  the  diverging  part  would  be  acceptable 
for  a certain  portion,  but  would  eventually  lead  to  a tie-up 
of  further  expansion  (within  the  constraints  of  the  formula- 
tions used,  like  continuous  expansion,  no  shock  waves  admitted). 

By  allowing  a change  in  the  cone  angle  6 along  the  axis,  which 
means  a departure  from  the^ straight  cone  geometry,  orderly  ex- 
pansion can  be  obtained.  The  choice  of  the  angle,  however,  is 
critical  since  too  small  an  angle  leads  to  poor  local  efficiencies 
(coupled  with  a drop  of  the  slip  ratio  K and  a decrease  of  the 
value  of  the  pressure  gradient  and  an  increase  in  droplet  size  d. 

Since  the  behavior  of  K,  p and  pressure  gradient  is  interrelated, 
it  is  possible  to  find  a suitable  nozzle  contour  (or  cone  angle 
6)  by  prescribing  a certain  gradient  dK/dT  at  the  inlet  of  each 
step.  A solution  obtained  for  a small  nozzle  is  shown  in 
Figure  10  where  dK/dT  was  uniformly  equal  to  0.02  at  the  begin- 
ning of  each  step. 

The  investigation  of  two-phase  flow  in  an  impulse  turbine  was 
necessarily  brief  and  is  given  in  Part  II.  The  results  given 
at  the  end  of  Part  II  show  that  the  relative  deviation  of  the 
droplets  from  the  steam  path  will  increase  with  the  Mach  index 
of  the  steam  flow,  the  ratio  of  droplet  diameter  to  the  radius 
of  curvature  of  the  flow  in  the  turbine  bucket  and  the  ratio  of 
droplet  diameter  to  the  "mean  free  path  length  of  the  steam  mol- 
ecules". Also,  the  larger  the  specific  volume  of  the  steam,  the 
larger  the  droplet  flow  deviation.  Combining  the  mean-free  path 
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length  with  the  other  steam  conditions  gives  for  the  "Stokes 
Regime" 


VCa  * TS 
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a a 


Vbnark 
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which  has  the  form  of  a Reynolds  Number,  with  d2/^  as  refer- 
ence length.  To  minimize  droplet  impingement  on  the  buckets, 

£.  large  radius  of  curvature  r^  of  the  flow  path  is  needed, 
together  with  low  steam  velocities  and  small  droplet  diameters. 

Even  with  droplet  impingement  and  liquid  film  formation, 
acceptable  performance  seems  possible  with  dished  radial  out- 
ward buckets  as  were  proposed  in  Reference  8. 

Further  studies  are  required,  especially  with  larger  machine 
dimensions.  Systematic  applications  of  the  computer  programs 
developed  should  give  nozzle  designs  of  the  characteristics 
required  for  high  performance  machines.  That  conclusion  is 
supported  by  the  results  of  Reference  7. 
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r~1-,  r-jv 


NOMENCLATURE 
(Page  1 of  4) 


Arabic 


Acceleration 

Cross-sectional  area  of  nozzle 
* Throat  area  of  nozzle 
*C"  Specific  heat  of  saturated  steam 

*C’  Specific  heat  of  saturated  water 

Index  for  independent  variable 
Velocity  of  steam 

Velocity  of  water  droplets 


Drag  coefficient 


v/2  (i  - i) 

: m0  m 

, Velocity  from  isentropic  expansion 

Droplet  diameter 
Nozzle  diameter 
, Droplet  drag 

= \fpv  Reference  velocity 

Index  for  successive  discretizations 
Static  enthalpy 

i Static  enthalpy  of  mixture 

Index  in  solution  matrix  T[J,I] 

Index  in  solution  matrix  T[J,I] 
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NOMENCLATURE 
(Page  2 of  4) 

K ^a^Cb  S^P  rati° 

L Power  loss 

m Mass  flow 

mT  Total  mass  flow 

N[H]  Number  by  which  step  size  is  divided  by 
p Static  pressure 

Q 4x (1  -x)  (n/k  - 1/v/K  ) 

rc  Latent  heat  of  evaporation 


Radius  of  curvature  of  path 
dv 

R Z dx/dp  + — — — 

Vv/K 

Re  Reynolds  number  dw/v  n_ 

3 3 

1 c 

s Entropy 

S 2dD/Ddp 


t Temperature  [°C] 
T Temperature  [°K] 
v Specific  volume 


NOMENCLATURE 
(Page  3 of  4) 


w 


Relative  velocity  between  droplet  and  steam 


W 

We 

/ 

X 

x* 

X 

X 


3 „ (1-x)  vb  (K- 1) * D 

'T  CD  i va  K i tans 

w2d 

2v  a 
a 

Vapor  mass  fraction  (quality) 
Axial  length  coordinate  of  nozzle 
xK  + (l-x)/K 
Xv/K  + (l-x)/v/K 


Y[C,H]  Dependent  variable,  approximation  of 

v. 


f-ui:  * (#,  ~ vb'gl 


Greek 

o Static  void  fraction 

6 Half  angle  of  conical  nozzle  element 

na  Viscosity  of  steam 

nN  Nozzle  efficiency  (local  and  overall) 

Xa  vana^e  "Mean  free  path  length  of  steam  molecules" 

p Mass  density 

a Surface  tension 
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NOMENCLATURE 
(Page  4 of  4) 


[ 


Subscripts  and  Superscripts 

0 Stagnation  conditions 

1 Nozzle  inlet 

a Steam 

b Liquid 

* Throat  conditions 


I 

i 

i 

I 


X 
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APPENDIX  II 


PROPERTIES  OF  WET  STEAM 
UP  TO  2 50°C 


1*  The  Callendar  Equation  of  State 


TKe  Callendar  equation  of  state  was  used 


va  = T ' °‘°75 


+ % ■ 


(ii-D 


where  T is  in  degrees  Kelvin,  and  the  specific  volume  in  m3/s. 


2.  The  Change  in  Vapor  Specific  Volume 


The  change  in  dv  is 

3. 


■ (tt)  dT  * (^)dP 


(II-2) 


From  Eq.  (Al)  the  partial  derivatives  are  as  folic 


(£) 


= £ + 0.075 


10  / 273 . 2 \ / 2 75 . 2 \ 

3 \~) 


(II-3) 
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(I 1-4) 


The  Clausius-Clapeyron  Eq . (60)  gives  the  required  dp/dT 


Ifr  - 


.V „ - V, 

a b 


(11-5) 


so  that  we  get  finally 


Va  _f3va\  dT  / 3va  \ 

Ip"  \Tr/p  3?  \-5p-/ 


CII-2') 


• fo)„  * ft)T  S ’ 


(II-2") 


j.  Saturation  Pressure,  Latent  Heat,  Enthalpy  and  Entropy  of  Liquid 

The  remainder  of  the  relations  required  was  compiled  from  Ref.  16 
and  with  exception  of  the  specific  heats,  is  summarized  in  Table 
1 1 - 1 - Single  primed  quantities  represent  saturated  water  proper- 
ties; double  primed  quantities  represent  saturated  steam  properties. 

4 . Specific  Heat  of  Saturated  Water  and  Steam 

The  specific  heat  of  the  saturated  liquid  is  calculated  using 
V.  Regnault  data,  see  Ref.  16,  page  227 


C.  = 4186(1  + 4x1 0 ' 5 1 + 9xl0~7t2) 


The  temperature  t is  to  be  taken  in  degrees  Celsius, 


(11*6) 


The  specific  heat  of  the  saturated  vapor  was  determined  from  the 
basic  equations  of  thermodynamics  as  follows 


c = c = C - T 
a P Po 
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( 1 1 * 7) 
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with  v = 

a p 


/ 7 7 3 7 \ n 

I ~ y — I + vb  according  to  Callendar,  so  that 

(S)  • 

P 


where 


C = 0.075  m3/kg,  n = 10/3, 


C = 1996.7  J/kgC, 
po 


the  low  pressure  ideal  gas  specific  heat  of  water  vapor 


We  finally  get 


C = 1996.7  + 1.0833 

a 


7 \ic  /'3 

(^)  f 


, (11-14) 


where  the  units  of  p are  [N/m2],  and  those  of  T in  degrees  KelvinCk], 
5.  Calculation  of  dx.  The  Change  in  Vapor  Mass  Fraction 


The  total  differential  dx  follows  from 


. /ax\  , (dx\ 

- W * H 


(II-1SJ 


Following  Ref.  (5)  we  get  for  a two-phase  mixture  like  wet  steam 


( —\ 

VP  Is 

(a). 


(9s/ 9p) x 
(9s/9x)p 


(11-16) 


= S"  - S'  * 


(11-17) 


where  s"  - s'  is  the  entropy  change  upon  evaporation.  Since 


s = xs"  + (l-x)s' 


(»), 


= x + C1-*)  %- 


(11-18) 


(11-19) 


The  differential  ds"  is 


(11-20) 


For  an  assumed  ideal  gas  model  for  the  saturated  vapor  of  specific 
heat  c , the  Eq.  of  Gibbs  is 


Tds  = c^dT  - vdp  » 


(11-21) 


and  therefore: 


(¥) 


(11-22) 


(^)o  ' ^ 


(11-23) 


Use  of  the  Clapeyron  Equation 


dp  = s”  - S’ 
^ ~ v"  - v' 


(11-24) 


gives  after  substitution  of  Eq . (11-24),  (11-22)  and  (11-23)  into 
Eq.  (11-20) 


ds"  _ v"  . p (v"  - v') 

^ ' ' T 


(11-25) 
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which  is  ready  for  substitution  into  Eq.  (11-19). 


r 

'• 


For  the  saturated  liquid  of  specific  heat  c’ 


ds'  = ds'  dT  _ c'  (v"  - v’) 

■ -ar  ap  - t (s„  _ sl) 


(11-26) 


which  is  to  be  substituted  into  Eq . (11-19)  also. 


Final  substitution  of  Eq.  (11-17)  and  (11-19)  into  Eq.  (11-16) 


gives 


3x  ) = x (vM - 

,^P/e  Tfc". 


f— > ♦ — _1_CP] 
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c. 1 (v"  - V*) 

T (s’  - s')2 


x(c'  - c ) 
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The  second  term  in  Eq.  (11-15)  may  be  written 


/ 9 x \ ds  dp 

\3s 3p  3T 

From  the  definition  of  efficiency,  n.  ■ di/vdp,  we  get  from  the 


Gibbs  Equation,  Tds  = di  - vdp , 


1 - n 


, -Tds 

N v dp 
m K 


-V!  * V 


(11-29) 


(11-30) 
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so  that  finally,  from  Eqs.  (11-15)  and  (11-17) 


dx  / d) 

vmC1  ‘ "N1 

TT  \T\ 

r/s 

T(v"  - v') 

(II-31) 


6 . Viscosity  of  Saturated  Steam 

Tl^e  final  formula  used  is  of  the  Sutherland  type,  as  follows 

_ / _ V»/2 


273  + C / T V 
'o  "T  +~C  [773/ 


(11-32) 


where 


nQ  = 8.6843  x 10’  [kg/ms] 

— [fc] 

T = Temperature  [°K] 

Sample  values  were  calculated  according  to  the  latest  international 
formula,  see  Reference  17;  they  are  tabulated  below  together  with 
results  from  the  above  formula. 

VISCOSITY  OF  SATURATED  STEAM 


t 

C°C] 

International 
Formula  1975 
n [kg/ms] 

Sutherland 
Formula 
n [kg/ms] 

100 

1 . 2 2 7760x10' 5 

1 . 2285x10* 5 

110 

1 . 262300x10’ 5 

1 . 2637x10* 5 

120 

1.297068x10' 5 

1 . 2987x10" 5 

130 

1 . 331996x10' 5 

1 . 3336x10 ' 5 

140 

1 . 367034x10' 5 

1 . 3684x1 0 ' 5 

150 

1 .402147x10* 5 

1 . 4030x10* 5 

J 

7 . Surface  Tension  of  Water 

The  "International  Association  for  the  Properties  of  Steam" 
released  in  September  1975  the  following  formula  for  the  sur- 
face tension  of  water  (that  is  valid  between  the  triple  point 
and  the  critical  point) 

a = 0.2358  (l  - 57775-5)  [l  * 0.625  (l  - ^47^)]  (II'33) 

J? 

. where  T is  to  be  given  in  [°K],  and  the  units  of  surface  tension  a 

are  [N/m],  see  Reference  17. 
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Table  1 1 - 1 . Saturated  Steam  Properties 


1 1 


l * 
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o 


log 


OH 


l1 ' 

1*btl(1 

- *y  ~ v (’ 

L ' 1 

3] 

T/Tk 

1 ♦ fTk  (1 

- T/Tk) 

= 2 


(II-8) 


p/p.  * 10  z Good  up  to  critical  point 


1996.7  T 


0 . 0 7 5 (n+ 1 ) j ' V'|p  + 1-9441xlOs J/kgK 


or 


(II-9) 


i * + r , which  is  more  accurate  than  Eq . ( 7 ) for  saturated  steam 
c 


346. 6^1^  - T [kJ/kg]  Latent  heat  of  evaporation 


(11-10) 


10.1796  log 


(ttt)  ■ 8- 


6229x10" 4 (T-273) 


+ 1.8836x10"® (T-273) 2 kJ/kg°K 


4 . 186 (t  +2 . 0x10  5t 2 ♦ 3 . 0x1 0"  7 1 3 ) [kJ/kg] 


-7.3 


(H-ll) 

(11-12) 


1 + v’p  = q + pxlO  5 [kJ/kg]  Enthalpy  of  Liquid  (11-13) 


— - 0.075 
P 


(273^2)' 


0/3 


+ v’ 


(II-l) 


Eq.  (II-8)  By  Smith,  Keyes  6 Gerry,  1.934  (Plank,  R.,  p.  112) 
(11-10)  By  M.  Jakob,  1935  (Plank,  p.  125) 

(11-12)  Regnault  Data  for  Spec.  Heat  of  Liquid*  see  Eq.  (II-6) 

(11-11)  By  integration  from  (11-12) 

(11-13)  ®y  integration  from  (11-12) 

(II-l)  Callendar  Eq.  of  State 

(II-9)  From  General  Eq.  of  Thermodynamics  and  Eq.  (II-l). 
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